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Foreword 

The  material  in  this  report  represents  a contri- 
bution to  Nucleation  II , Marcel  Dekker , N.Y.  (in  press), 
edited  by  A.C.  Zettlemoyer.  This  contribution  was 
solicited  by  the  editor  who  also  assembled  Nucleation , 
Marcel  Dekker,  N.Y.,  1969  (now  called"Nucleation  I"). 

Since  the  material  includes  original  contributions, 
it  is  first  published  as  an  ONR  Report. 


ABSTRACT 


The  fundamentals  of  gasdynamics  of  one-dimensional 
flows  with  condensation  are  reviewed  paying  special  at- 
tention to  concepts  from  nonequilibrium  flows , homoge- 
neous nucleation,  droplet  growth.  Flow  systems  for 
experimental  research  in  condensation,  namely  steady- 
flow  Laval  nozzles,  Ludwieg  tubes,  shock  tubes  and  high 
intensity  supersonic  molecular  beams  are  discussed. 

New  and  previously  published  experimental  results  on 
homogeneous  nucleation  of  eight  vapors  , obtained  using 
gasdynamic  methods,  were  critically  examined  and  com- 
pared on  a uniform  basis  with  steady  state  homogeneous 
nucleation  theory.  It  is  found  that  experiments  with 
water,  nitrogen  and  ethanol  vapors  agree  with  the  classi- 
cal theory  of  homogeneous  nucleation  or  the  statistical 
mechanics  theory  based  on  the  models  of  Kuhrt  or  Dunning. 
This  study  also  revealed  that  in  a number  of  experiments 
with  several  vapors , in  which  much  higher  nucleation 
rates  than  those  calculated  from  the  classical  theory 
were  thought  to  have  been  observed,  involved  ambiguous 
interpretations  and  assumptions.  To  date,  no  experiment 
using  any  method  has  shown  conclusive  agreement  with  the 
Lothe  and  Pound  model  of  the  statistical  mechanical 
theory  of  homogeneous  nucleation.  Finally,  a bibliography 
of  condensation  studies  in  compressible  flow  systems  is 
provided  in  the  appendix. 
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1 HISTORICAL  INTRODUCTION. 

The  first  international  congress  wholly  devoted 
to  high  speed  flight  and  gasdynamics*  was  held  in  Rome 
in  1935.  At  that  meeting,  Prandtl  (1935)  showed  a 
schlieren  picture  of  nozzle  flow  (Fig. 1.1)  and 
he  pointed  to  the  then-mysterious  disturbances  seen  in 
the  supersonic  flow  to  the  right  of  the  nozzle  throat. 
Wieselsberger**  in  the  ensuing  discussion  suggested 
that  condensation  of  water  vapor  may  have  been  responsible 
for  these  shock-like  phenomena  since  Prandtl' s supersonic 
wind  tunnel  was  operated  with  moist  air.  Shortly  after 
the  Rome  meeting  this  guess  was  indeed  confirmed  at 
Wieselsberger ' s laboratory  in  Aachen  with  the  results 
published  later  by  Hermann  (1942).  Subsequently 

air  driers  were  introduced  in  supersonic  wind  tunnel 
circuits  to  avoid  condensation  disturbances  which  rendered 
quantitative  model  tests  impossible.  Thus,  condensation 
problems  entered  gasdynamics  and  aeronautics  by  an  acci- 
dental  discovery  while  water  condensation  has  long  ago 
been  known  to  occur  in  steam  nozzles.  These 

* Gasdynamics  is  the  regime  of  fluid  mechanics  where  the 
compressibility  of  the  fluid  needs  to  be  taken  into  account. 
**  For  a brief  history  of  these  early  events  and  detailed 
references  see  the  Appendix  of  Wegener  (1975). 
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2. 


convergent-divergent  nozzles  introduced  by  the  Swedish 
engineer  de  Laval  in  1883  (e.g.  Stodola  1927)  produced 
supersonic  steam  flow  in  steam  turbines.  The  appear- 
ance of  condensation  in  such  nozzles,  however,  was  de- 
layed with  respect  to  the  coexistence  line  as  shown  in 
Fig.  1.2.  From  many  such  observations  an  empirical 
condensation  curve  emerged  which  was  named  the  "Wilson 
line"  to  honor  the  early  investigator  of  delayed  con- 
densation in  cloud  chambers  (e.g.  Wilson  1897). 

Condensation  of  water  vapor  in  air  in  supersonic 
nozzles  obeyed  empirical  relationships  and  a function 

AT  , = f($  , dT/dt)  , (1.1) 

ad  o’  ’ 


could  be  established  experimentally.*  The  adiabatic 
supercooling  ATad  is  taken  to  be  the  temperature  differ- 
ence between  the  location  of  saturation  and  the  onset  of 
condensation  on  an  adiabatic  (here,  an  isentropic) 
expansion. 


* For  symbols  note  Section  8. 


- 
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This  term  is  also  illustrated  in  Fig.  1.2.  It  is  more- 
over apparent,  that  the  "critical  supersaturation" , i.e. 
the  ratio,  Pvjc/poo > of  the  partial  pressure  of  the  vapor  to  its 
saturation  pressure  at  the  onset  of  condensation 
is  an  equivalent  method  to  describe  the  state  of  the 
spontaneous  condensation  in  the  supersaturated  state. 

For  ATacj>0,  we  have  pv^/poo>l.  In  supersonic  flow  we 
normally  experience  high  values  of  the  critical  super- 
saturation, therefore  it  is  usually  easier  and  more 
accurate  to  express  the  onset  of  condensation  in  terms 
of  AT  ^ since  slight  errors  in  temperature  result  in 
large  errors  in  the  tabulated  saturation  pressure. 

The  primary  effect  on  the  value  of  the  supercooling 
in  Eq.(l.l)  is  due  to  the  initial  partial  pressure  of 
water  vapor,  or  the  relative  humidity  in  the  supply  of 
the  nozzle,  4>  . Experimental  data  on  this  effect  in  noz- 
zle flow  are  shown  in  Fig.  1.3  at  a cooling  rate  of 
about  10 5C/sec . Next,  Fig.  1.4  shows  the 

effect  of  a variable  cooling  rate  observed  in  a Prandtl- 
Meyer  expansion*  with  a functional  relationship  also  shown 

* Such  an  isentropic  expansion  arises  in  a steady  super- 
sonic flow  around  a sharp  convex  corner.  It  is  the  two- 
dimensional  analog  of  the  unsteady  expansion  in  shock 
tubes  to  be  discussed  in  Section  2c. 


4. 

to  admit  the  proper  physical  limit  at  ATac)  = C for  a 
vanishing  cooling  rate.  We  note,  for  this  experiment 
as  for  other  gasdynamic  methods,  that  the  cooling  rates 
are  appreciably  higher  than  those  of  other  techniques. 

It  took  until  1940  before  this  empirical  result 
and  in  fact  the  entire  gasdynamic  process  were  physi- 
cally understood.  Displaying  remarkable  insight, 
Oswatitsch  (1941,  1942)  provided  experimental  and  theo- 
retical evidence  on  the  nature  of  the  condensation 
process  in  supersonic  flow  (homogeneous  nucleation, 
growth  laws,  time  scales,  heat  addition,  etc.)  that 
still  stands  today.  Moreover,  he  was  the  first  to 
suggest  supersonic  flow  methods  to  study  condensation 
itself  by  proposing  to  measure  the  surface  tension  of 
small  clusters  indirectly.  This  task  still  eludes  us 
today.  Thus,  he  added  a new  tool  to  the  chemical 
physics  laboratory  that  rivals  Wilson's  cloud  chamber 
in  importance  (Dunning  1960). 
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Finally,  Oswatitsch  noted  the  cross  connection  of  con- 
densation in  supersonic  wind  tunnels  with  that  in  steam 
turbine  nozzles  and  he  employed  identical  theoretical 
methods  for  the  understanding  of  both. 

Gasdynamic  methods  have  since  been  enlarged  by 
adding  shock  tubes  and  other  devices.  However,  they 
all  have  some  common  features  which  set  them  apart  from 
other  techniques.  For  one,  heterogeneous  nucleation  on 
aerosols  present  in  the  flow  is  unimportant.  The  short 
time  scales  of  supersonic  flow  are  such  that  even  maxi- 
mized rates  of  condensation  on  foreign  nuclei  are  negli- 
gible in  relation  to  condensation  produced  by  homogeneous 
nucleation.  This  situation  is  unchanged  even  in  the 

i 

presence  of  large  numbers  of  aerosols,  say  106/cm3, 

(Oswatitsch  1942,  Wegener  1969).  It  is,  in  fact,  dif- 
ficult to  affect  the  normally  homogeneous  nucleation  by 
seeding  the  moist  air  in  the  supply  of  a supersonic  noz- 
zle (Buckle  and  Pouring  1965).  Supersaturation  is  gen- 
erally found  to  be  high  owing  to  the  high  cooling  rates 
because  the  short  time  scale  permits  only  a reduced 
number  of  collisions  of  the  condensing  vapor  in  relation, 
for  example  to  a cloud  chamber.  However,  these  effects 
still  permit  a detailed  resolution 
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of  the  condensation  process.  Rather  than  the  single 
measurement  of  the  critical  supersaturation  possible 
in  Wilson  and  diffusion  cloud  chambers,  condensation 
rate,  nucleation  rates,  droplet  growth,  etc.,  may  be 
measured  quantitatively  in  a precisely  prepared  en- 
vironment . 

2 COMPRESSIBLE,  INVISCID,  ONE-DIMENSIONAL  FLOW  OF 
PERFECT  GASES.* 

a)  Basic  Considerations 

Effects  of  compressibility  arise  in  fluid  mechanics 
if  flow  speeds  exceed  about  one  third  of  the 
speed  of  sound,  and  this  technologically  important 
regime  has  received  much  attention  in  the  last  forty 
years  or  so.  The  advent  of  supersonic  flight,  rocketry, 
atmospheric  entry  from  space,  etc.  has  moreover  spurred 
much  recent  progress.  In  such  flows, the  density  varies 
and  in  addition  to  the  equations  of  motion, an  equation 

* This  section  is  addressed  to  those  who  may 
not  be  familiar  with  gasdynamics.  For  details  we  refer 
to  textbooks  such  as  those  by  Liepmann,  ri.W.  and  Roshko, 
A.,  Elements  of  Gasdynamics , Wiley,  New  York  1957, 
Becker,  E.  , Gas  Dynamics , Academic  Press,  New  York  1968  , 
and  Vincenti,  W.G.  and  Kruger,  C.H.,  Introduction  to 
Physical  Gasdynamics , Wiley,  New  York  1966. 
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of  state,  p = p(p,  T),  is  required.  Conversely,  gas- 
dynamics  had  a strong  impact  on  chemistry  (e.g.  Wegener 
1966)  since  well  understood  gasdynamic  processes  such 
as  shock  waves  provide  tools  for  the  study  of  reactions 
at  controlled  time  resolutions  that  have  traditionally 
not  been  available.  Nonequilibrium  flows  are  important 
(e.g.  Becker  1972)  where  one  or  more  modes  of  molecular 
motion,  etc.  of  a gas  mixture  depart  from  thermodynamic 
equd librium.  Rapid  condensation  by  homogeneous  nucle- 
ation  represents  such  a nonequilibrium  situation. 
However,  before  studying  real  gas  flows,  a few  simple 
gasdynamic  concepts  need  to  be  clarified.  This  dis- 
cussion will  set  the  stage  for  an  understanding  of  the 
particular  advantages  of  high  speed  flow  experiments  in 
relation  to  the  more  traditional  devices  in  the  field. 
Gasdynamics  is  described  in  relation  to  condensation  in 
greater  detail  by  Stever  (1958),  Wegener  and  Mack  (1958) 
and  Wegener  (1969). 
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it  is  found  in  practice  that  at  low  shear  rates , the 
internal  friction  (or  turbulent  motion)  can  often  be 
neglected  and  the  fluid  can  be  treated  as  if  it  were 
inviscid.  This  idealization  of  the  real  fluid  is 
never  valid  near  solid  walls.  In  a narrow 
layer,  the  so-called  boundary  layer,  the  flow  speed  is 
reduced  from  the  undisturbed  free  stream  speed  to  zero 
velocity  at  the  wall.  The  structure  of  this  boundary 
layer  (laminar  or  turbulent) , and  its  lateral  extent  are 
governed  by  the  well-known  Reynolds  number.  This  di- 
mensionless similarity  parameter  is  given  by 


D _ pui,  uJl 

Re  = V = ~ > 


(2.1) 


expressing  the  ratio  of  the  inertial  to  the  viscous 
forces  acting  on  a fluid  element.  Generally  at  Reynolds 
numbers  of  the  order  of  10**  to  10s  as  found  in  most  con- 
densation experiments  in  gasdynamics , the  boundary  layers 
near  the  walls  of  a nozzle  are  turbulent. 

Next,  we  recall  the  definitions  of  the  Mach  number: 

M = u/a  , (2.2) 

and  the  Knudsen  number,  relating  the  mean  free  path  to 
a characteristic  length: 


Kn  = A/4. 


(2.3) 
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By  combining  Eqs.  (2.1),  (2.2)  and  (2.3)  and  making 
use  of  concepts  from  kinetic  theory  we  find 


where  y is  the  ratio  of  the  specific  heats.*  Here  all 
three  important  similarity  parameters  of  compressible 
flow  appear.  Excepting  high  intensity  molecular  beams 
produced  by  supersonic  nozzle  flow,  large  Knudsen 
numbers  characteristic  of  the  free  molecule  flow  (FMF)  regime 
are  not  found  in  gasdynamics  research  on  conden- 

sation. Rather, Kn<<l  for  the  previously  quoted  high 
Reynolds  numbers,  and  continuum  flow  prevails.  At  high 
Reynolds  numbers , viscous  effects  are  small  and  the 
boundary  layers  are  found  to  be  very  narrow  (e.g.  Fig. 

8,  Wegener,  1969).  For  such  nearly  inviscid  flows5 a 
small  empirical  correction  can  be  applied  to  the  geo- 
metrical nozzle  cross-sectional  area  resulting  in  an 
"effective  area"  to  be  used  in  the  inviscid  flow  treat- 
ment to  be  discussed  next. 

One-dimensional  Flow:  The  requirement  of  low  shear 

* It  is  customary  in  gasdynamics  to  use  the  symbol,  y, 
for  this  term.  We  shall  therefore  later  denote  the 
surface  tension  by  a , in  contrast  to  other  chapters  in 
this  volume. 
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to  produce  practically  inviscid  flows  further  demands 
flow  systems  of  low  or  no  wall  curvature  in  addition 
to  the  previous  criterion  of  high  Reynolds  numbers. 
Typical  nozzles  used  for  condensation  studies,  again 
with  the  marked  exception  of  gasdynamic  molecular  beams 
(e.g.  Anderson  et  al . 1966),  consist  of  converging- 
diverging  ducts  having  flat  walls  with  total  included 
angles  of  about  5 to  10  degrees.  These  plane  walls 
are  joined  by  cubic  curves  in  the  throat  region  to 
provide  a continuous  second  derivative  of  the  surface 
profile.  Consequently,  a uniform  acceleration  of  the 
flow  at  the  junctures  of  the  cubic  curve  and  the  plane 
is  found.  In  such  ducts  the  concept  of  one-dimensional 
flow  is  applicable  with  a practically  uniform  velocity 
distribution  across  the  channel,  i.e.  at  a right  angle 
to  the  streamlines.  In  fact  if  we  have  A(x)  given  by 
a static  pressure  calibration  of  a nozzle  (see  Section 
2b),  we  can  later  find  u(x)  from  the  equation  of  motion. 

In  turn  in  shock  tubes  of  constant  cross  section 
the  flow  is  one-dimensional  by  definition  provided  again 
that  the  Reynolds  number  is  high  enough  to  insure  a 


11 


negligible  boundary  layer  thickness.  In  sum,  the 
treatment  of  nozzle  and  shock  tube  flows  at  high 
Reynolds  numbers  and  low  supersonic  Mach  numbers  will 
permit  us  to  apply  the  equations  of  motion  simplified 
for  inviscid  (ideal)  and  one-dimensional  flow  res- 
pectively . 

Equation  of  State:  Finally,  a substantial 

analytical  advantage  is  obtained  if  the  gases , vapors 
and  their  mixtures  can  be  treated  as  thermally  and 
calorically  perfect  gases  until  condensation  appears. 

This  assumption  applies  in  particular  to  gasdynamic 
experiments  designed  to  study  the  condensation  process 
itself.  Exceptions  are  found  for  thermal  non-ideality 
at  higher  pressures  in  steam  nozzles  (e.g.  Barschdorff 
1971).  The  need  for  real  gas  equations  including  those 
for  the  supersaturated  states  becomes  particularly  acute 
in  modern  high  pressure  steam  turbines  since  states  at 
extreme  conditions  of  say  50  bar  pressure  at  condensa- 
tion are  required  to  make  efficient  use  of  nuclear  plants 
(Gyarmathy  et  al.  1973).  On  the  other  hand,  in  shock 
tube  work,  the  study  of  condensation  of  metal  vapors 
(e.g.  Kung  and  Bauer  1971) 


may  require  that  caloric  non-idealities  appearing  at 
high  temperatures  are  included. 

However,  aside  from  these  less  frequent  situations, 
the  ideal  gas  assumption  is  valid  and  the  equation  of 


state  is  given  by 


p = RT  E ( P • /y  • ) , 

j 3 3 


(2.5) 


where  the  subscript  j denotes  the  j 'th  species.  Eq. 
(2.5)  can  also  be  written  as 


p = p(R/y)T  , 


(2.6) 


where 


P = ^ , 


(2.7) 


1/W  = E (a)j/yj)  , 


(2.8) 


using  the  mass  fraction,  uk  = Pj/p. 
of  the  specific  heats,  y,  is  given  by 


The  ratio 


c Eu>  .c  . 

= _£  = 3_£2 

c Euj.c  . * 

v J V] 


(2.9) 


It  is  practical  in  gasdynamics  to  express  heat  capacities, 
enthalpy,  etc.  in  terms  of  unit  mass  and  to  employ  lower 
case  letter  symbols  to  designate  this  fact.  Finally, 


Eqs . (2.5)  and  (2.6)  can  later  readily  be  adapted  to 
include  partially  condensed  flow  by  introducing  the  mass 
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fraction  of  the  condensate,  g = p^/p^.  Here,  p^  is  the 
density  of  condensate  referred  to  the  volume  of  the  gas 
phase  and  p^  is  the  total  density  of  the  two-phase  med- 
ium (see  Section  5a). 

For  a thermally  and  calorically  perfect  gas,  the 
acoustic  velocity,  defined  by: 

a = /Op/3p)s  , 
is  given  by 

a = /y (R/p )T  , 
with  y and  y given  by  Eqs. 

Using  Eqs.  (2.11)  and  (2.2), 
ideal  gases : 

M = — • (2.12) 

✓y(R/p)T 

b)  Nozzle  Flow. 

Isentropic,  one-dimensional,  steady*  nozzle  flow 
exemplifies  the  special  attractiveness  of  gasdynamic 
techniques  in  condensation  studies.  With  the  assumptions 
of  the 


(2.10) 

(2.11) 

(2.9)  and  (2.8)  respectively, 
we  have  the  Mach  number  for 


* Steady  flows  are  those  in  which  no  parameter  changes 
with  time  at  a given  location. 
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last  section,  the  equations  to  be  derived  for  this  flow 


willf  moreover^  be 


useful  for  a determination 


of  the  starting  point,  or  "onset",  of  condensation  in 
nozzleSj since  up  to  this  point  the  condensation  pro- 
cess cannot  yet  interfere  with  a supersonic  flow.  The 
results  are  as  follows.* 

The  conservation  of  mass  is  expressed  by : 

puA  = Q = constant,  (2.] 


d_p  + du  + dA  _ q 
P u A 


Newton's  second  law  of  motion  applied  to  a moving, 
compressible  fluid  yields  Euler's  equation 

u du  + ^ = 0 . ( 2 . 

P 

Similarly, for  later  use  in  one-dimensional  flow  in  the 
x-direction,  Eqs.  (2.13)  and  (2.15)  can  be  written  as 

d(PuA>  =0  (2  i 

dx  U ’ 


du  , 1 dp  n 

u -7—  + — = 0 , 

dx  p dx 


* For  details  note  Wegener  and  Mack  (1958),  and  Wegener 
(1969)  in  addition  to  the  textbooks  cited. 


15 


where  the  term  u(du/dx)  expresses  the  convective  acceler- 
ation in  steady  flow.  We  realize  that  while  in  steady 
flow  no  flow  parameter  changes  with  time,  a gas  sample 
undergoing  the  expansion  alters  its  properties  from  point 
to  point  since  it  is  convected  to  different  locations 
having  different  states. 

By  combing  Eqs.  (2. 6), (2. 12), (2. 14),  and  (2.15)  we  can 
derive  the  well-known  area-velocity  relation  for  compress- 
ible flow, 


du  -dA/A 
u " 1 - M2  ‘ 


(2.18) 


The  fundamental  equation  for  incompressible  flow  yields 
the  well-known  result  du/u  = - dA/A  since  M = 0.  However, 
if  M i 0,  the  following  pattern  emerges.  For  0 < M < 1 
(or  subsonic  speeds  with  u<a),  Eq.  (2.18)  behaves  qualita- 
tively as  in  incompressible  flow;  a decrease  in  area  re- 
sults in  an  increase  in  velocity.  For  M > 1 (or  super- 
sonic speeds  with  u > a) , the  denominator  of  Eq.  (2.18) 
becomes  negative,  resulting  in  a positive  coefficient  for 
dA  on  the  right hand  side  of  the  equation.  Now  an  increase 


4 


1 
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of  flow  speed  is  caused  by  an  increase  in  area,  in  con- 
trast to  the  low  speed  situation.  Finally  at  M = 1 , 

(or  at  sonic  speed  with  u = a),  we  have  (1-M2)  = 0 in 
Eq.  (2.18).  In  steady  flow  this  condition  can  only  be 
realized  at  a "throat"  where  dA/A  = 0,  because  du/u 
must  remain  finite.  In  sum,  flow  from  a low  speed  re- 
servoir through  a converging-diverging  nozzle  can  attain 
sonic  speed  at  the  throat  and  subsequently  proceed  to 
supersonic  speeds  in  the  diverging  passage. 

To  obtain  the  local  properties  of  the  nozzle  flow, 
we  add  the  energy  equation  for  steady  adiabatic  flow: 

dh  + udu  =0.  (2.19) 

With  the  thermodynamic  assumptions  of  the  previous  section 
we  find  for  the  enthalphy 

h = c T + const.  , (2.20) 

P 

and  Eq.  (2.19)  can  be  integrated  to  give 

c T + u2 / 2 = c T , (2.21) 

P P o ' 

where  the  subscript,  o,  refers  to  the  known  nozzle  supply 
or  reservoir  state  where  u = 0.  Moreover,  Tq  may  be  de- 
fined as  the  stagnation  temperature  of  the  gas  after  it 


* 
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has  been  brought  to  rest  adiabatically . With  the  pre- 
vious assumptions  and  definitions  we  find: 

T_  = (j.  ♦ 1 . <2 


In  an  isentropic  expansion: 


18. 


t-w 


1/ (y-1) 


(2.28) 


Equation  (2.27)  tells  us  that  to  achieve  sonic  speed  at 
the  throat,  a pressure  ratio  across  the  nozzle  of  at  least 
pQ/p*  = 1*9  is  needed  for  y = 7/5  = 1.4.  Indeed,  Eq.  (2.18) 
does  not  prohibit  subsonic  flow  at  the  throat  at  insuffi- 
cient pressure  ratios  across  the  nozzle  with  a subsequent 
deceleration  in  the  diverging  part  of  the  nozzle.  Moreover, 
should  the  pressure  ratio  exceed  that  given  by  Eq.  (2.27), 
the  constant  in  Eq . (2.13)  can  now  be  determined  by: 


* * * * 

puA  = p a A = Q = const  . , 

ft  V 

from  a fixed  supply  since  a /a  = (T  /T  ) , 


we  find 


* * 

A _ p a 
A^  ~ pu  ’ 


(2.29) 


From  Eq . ( 2 . 29 ) 


(2.30) 


and  from  the  formulas  derived  so  far,  the  nozzle  area  ratio 
as  a function  of  Mach  number. 


£ - «[(¥)/( 


+ 


2 (y-1 i 


(2.31) 


may  be  obtained.  The  function  of  Eq.  (2.31)  is  double- 
valued and  for  each  area  ratio  there  are  two 


r 


4i  i .uunum 
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values  of  the  Mach  number  for  subsonic  and  supersonic 
speed  respectively,  while  A/A*  = 1 at  M = 1 . 

These  results  can  now  be  put  to  practical 
use.  For  a given  converging-diverging  nozzle  with  known 
geometry,  we  have  A(x)  and  thus  A/A*(x).  From  Eq.  (2.31)5 
we  next  find  M(x)  and  the  previous  formulas  yield  T/T  , 
p/pQ,  and  p/pQ  at  each  location.  A single  measurement 
in  a supersonic  nozzle  flow  at  a known  location  there- 
fore fully  defines  all  flow  parameters.  These  parameters 
can  be  tabulated  as  a function  of  Mach  number  as  given 
for  y = 7/5  =1.4  (diatomic  gases)  in  Appendix  1,  and 
they  are  shown  in  Fig.  2.1.  A rapid  drop  of  pressure 
with  Mach  number  is  apparent  and  relatively 

high  values  of  the  supply  pressure  are  required  for  the  flow 
to  remain  in  the  continuum  regime  at  higher  Mach  numbers. 

For  water  vapor,  ethanol  and 

many  other  vapors,  however,  saturation  is  usually  attained 
at  M < 1,  ahead  of  the  throat.  The  delayed  condensation 
occurs  usually  at  low  supersonic  Mach  numbers  (M  < 1.5),  or 
shortly  past  the  throat.  Therefore  even  at  relatively 
low  supply  pressures,  free  molecule  flow  does  not  yet 
enter.  The  geometry  of  the  nozzle  throat  region 


controls  the  expansion  rate  and  thus  the  cooling  rate 
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of  the  vapor  sample.  In  steady  nozzle  flow  such  time 
dependent  variables  can  be  found  from  the  local  grad- 
ients by 


d(  ) _ „ d(  ) 
dt  ~ dx  * 


(2.32) 


where  any  of  the  derived  parameters  may  be  put  in  the  paren- 
theses. For  a two-dimensional  nozzle,  the  change  of 
area  ratio  with  distance  in  the  throat  region  can  be  re- 
lated geometrically  to  the  radius  of  curvature  of  the 
throat,  R* , and  the  nozzle  throat  height,  h*  (Wegener 
1964).  We  then  find 


d (T/T  ) 
o_ 

dt 


it  it 

(R  h ) 


Y-l 


[ ( Y + l ) / 2 ] ' 


(2.33) 


or 


j~d(T/To)  f _ 


it  it 

Ca  (R  h ) 
o 


(2.34) 


For  diatomic  gases,  y = 7/5,  and  C = 0.278.  From  Eq. 
(2.34)  we  find  remarkably  high  cooling  rates  of  the 
order  of  10 6 C/sec  for  small  research  nozzles  (See 
Section  6). 


c)  Shock  Tube  Flow 

A simple  shock  tube  is  a duct  of  uniform  cross  sec- 
tion in  which  a removable  diaphragm  divides  two  chambers. 
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In  these  chambers  there  are  carefully  prepared  gases  or 
gas  mixtures  at  different  pressures  and  temperatures 
prior  to  a given  experiment.  These  gases  undergo  con- 
trolled changes  in  a wave  system  resulting  from  the 
rupture  of  the  diaphragm  to  equalize  conditions  in  the 
two  chambers  as  shown  schematically 

in  Fig.  2.2.  Condensation  research  has  utilized  several 
different  sections  of  the  wave  system  to  be  des- 

cribed in  Section  4b.  The  gasdynamics  of  shock  tube  flows 
is  discussed  in  detail  in  the  textbooks  cited  and 
special  applications  to  physical  chemistry  exploiting 
the  well  controlled  changes  in  the  gas  sample  are  given 
e . g . by  Bradley . * 

After  diaphragm  rupture  a shock  wave  develops 
rapidly  and  this  wave  proceeds  at  a uniform  speed  into 
the  undisturbed  low  pressure  gas  at  condition  1.  The  shock  Mach 
number  is  defined  by  M = u /a,  where  the  subscript,  s, 

S o _L 

relates  to  the  speed  of  the  shock  wave.  We  have  from  Eq. 

(2.11),  a^  = /y^ (R/y  ^ )T^ , and  the  pressure  jump  across 
the  wave  is  given  by  the  well-known  Rankine-Hugoniot 
relation  for  a calorically  and  thermally  perfect  gas  by; 

pi  ■ 1 * ^ • (2-35) 


* John  N.  Bradley,  "Shock  Waves  in  Chemistry  and  Physics," 
Methuen  Co.,  London,  Wiley,  N.Y.  1962  . 
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where  the  subscript  2 denotes  the  state  at  higher  pres- 
sure and  temperature  behind  the  shock.  As  seen  in  Fig. 
2.2,  Pi+>P^5  and  while  the  shock  moves  to  the  right,  an 
isentropic,  unsteady  expansion  fan  proceeds  into  the 
high  pressure  chamber.  The  front  of  this  expansion  moves 
at  the  sound  speed  of  the  medium  ahead: 


a4= /Y4(R/y4)T4  , 


(2.36) 


where  again  as  in  all  following  equations , y and  y are 
given  by  Eqs.  (2.8)  and  (2.9).  The  expansion  fan  shown 
in  Fig.  2.2  is  ideally  centered  at  the  diaphragm  location 
at  time  zero,  and  each  straight  line  in  the  fan,  being  a 
characteristic  of  the  governing  partial  differential 
equation,  defines  wave  propagation  at  constant  speed,  c, 
by 

c = x/t  = u - a.  (2.37) 

Here,  u,  is  the  local  flow  speed  of  the  gas  moving  towards 
the  lower  pressure.  At  the  head  of  the  expansion  fan, 


u = 0,  and  c = a^.  At  the  tail  of  the  same  fan  where 


c3  = x/t  = u3  - a3  ’ 


(2.38) 


the  outflow  speed  of  the  gas  initially  located  in  the  high 
pressure  chamber  may  be  locally  subsonic,  sonic  or  super- 
sonic depending  on  the  slope  (x/t)  of  the  last  character- 


istic Cg , i . e . 


a3  < U3  ’ 


(2.39) 


23 


It  can  be  shown  that  the  propagation  speed  of  this  last 
characteristic,  i.e.,  the  speed  of  the  tail  of  the  fan 
is  given  by: 


The  pressure  on  each  characteristic  in  the  expansion 
fan  may  be  related  to  the  corresponding  temperature 
and  density  in  the  isentropic  expansion  by  Eq. 

(2.23). 

The  local  wave  speed  can  be  expressed  in  terms  of  the 
local  flow  speed  and  the  conditions  in  region  4 by: 

Y4  + 1 

c = - a4  + — 2 u . (2.41) 


This  equation  allows  us  to  determine  the  path  of  an  actual 
particle  (or  small  volume)  of  the  flow  starting  at  some 
initial  location  as  shown  in  Fig.  2.2.  This  particle 
path  x = x(t)  is  obtained  by  integrating  the  local  fluid 
velocity  in  Eq.  (2.41): 


dx  _ 

at  = u 


y4+i  t 


^ + 


Yii+1  a4 


(2.42) 


The  history  of  a particle  can  now  be  traced  in  the  x - t 
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r 


it 


plane  with  initial  particle  position,  xq  at  time 

t (t  = x / a,  ) , by: 
o o o 4 J 


x(t) 


(2.43) 


Returning  to  Fig.  2.2  we  note  a last  feature  of  the 
flow,  the  so-called  contact  surface  separating  states  2 
and  3 corresponding  to  the  gases  initially  in  either 
chamber.  The  contact  surface  may  also  be  viewed  as  a 
virtual  piston  set  abruptly  into  uniform  motion.  We  find 
u2  = u3  anc*  ^2  = P3  * ^>u"t  ^3  < ^2  s^-nce  ‘t^ie  expanded  gas  on  the 
left  has  been  cooled  while  the  shocked  gas  on  the  right  has 
been  heated.  The  contact  surface  thus  represents  an 
entropy  discontinuity. 

Combining  the  conservation  laws,  we  finally  arrive 
at  the  shock  tube  operating  equation 


£4 

Pi 


Pi 


1 - 


( Yu-1)  (a-i /au) 


(p0/pn ) 


A 


t/2Yi+(Y1+1^P2/Pl)"1" 


-(2y4)/(y4-1) 

, (2.44) 


where  we  can  insert  Eq.  (2.35)  for  p2/p^-  This  substi- 
tution links  the  shock  Mach  number  directly  to  the  ini- 
tial conditions  that  are  known  a priori  and  thus  the 
remaining  parameters  can  be  computed.  The  maximum 


expansion  and  correspondingly  the  lowest  pressure  and 
temperature  respectively  achieved  in  the  isentropic, 
unsteady  fan  is  found  e.g.,  from  P3/P4  = (P2^Pi ^ (P4/P-L ^ 
in  conjunction  with  Eq.  (2.44)  since  p^  = pj.  Instanta- 
neous distributions  of  these  properties  at  time  t^,  are 
shown  schematically  in  the  lower  part  of  Fig.  2.2  to- 
gether with  symbols  denoting  the  conditions  in  the  tube. 

At  high  Reynolds  numbers,  small  viscous  effects  or 
the  minor  consequences  of  finite  diaphragm  rupture  can 
be  accounted  for  as  will  be  seen  in  Section  4.  The 
reasons  for  this  fortunate  state  of  affairs  were  outlined 
in  Section  2a.  However,  a check  of  the  ideality  of  the 
tube  operation  is  usually  provided  by  the  redundancy 
inherent  in  the  measured  variables  entering  the  equations 
of  motion.  Usually  the  shock  speed,  ug , is  measured  and 
thus  P2/P1  can  computed  directly  from  Eqs.  (2.12)  and 
(2.35).  Independently, with  known  initial  conditions, 

Pg  = Pj  is  measured  and  overall  agreement  with  Eq.  (2.44) 
can  thus  be  checked.  For  the  application  of  these  flows 
to  condensation  studies,  we  refer  to  Section  4. 

3 REAL  GASES  AND  FLOW  WITH  CONDENSATION 
a)  Nonequilibrium  Flow 

In  certain  flows , the  behavior  of  real  gases  may  de- 
viate from  that  of  ideal  gases.  In  terms  o 
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the  equation  of  state,  such  deviations  may  generally  be 
divided  into  two  classes:  those  resulting  in  an  equation 
of  state  involving  two  independent  variables,  e.g.,  for 
the  internal  energy, 

e = e (p  , s) , (3.1) 

and  those  requiring  more  than  two  independent  variables, 
such  as 

e = e(p,  s,4>).  (3.2) 

Here  the  additional  so-called  progress  variable,  <J> , is 
required  to  fix  the  thermodynamic  state  of  the  gas.  Non- 
idealities of  the  kind  of  Eq.  (3.1)  have  been  discussed 
before  in  conjunction  with  high  pressure  steam.  No  con- 
ceptual difficulty  arises  in  the  gasdynamic  treatment 
of  the  flow  since  the  function,  Eq.  (3.1),  may  be  given 
in  the  form  of  diagrams,  tables,  or  simply  in  more  compli- 
cated equations  than  the  ones  for  the  ideal  gas. 

Equation  (3.2)  on  the  other  hand  applies  to  situations 
where  departures  from  the  thermodynamic  equilibrium  appear. 
This  may  happen  if  chemical  reactions,  excitation  or  de- 
excitation of  some  internal  degrees  of  freedom  of  molecules, 
dissociation,  ionization  or  changes  of  phase  occur.  The 
progress  variable,  $ , must  therefore  be  defined  differently 
in  each  specific  non-equilibrium  situation.  When  the  flow 
is  in  the  continuum  regime,  i.e.,  when  Kn<<l,  the  gas  mole- 
cules may  be  considered  to  be  in 
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translational  equilibrium,  and  the  pressure  and  a "trans- 
lational temperature"  may  be  defined  for  the  gas  without 
ambiguity.  Moreover,  at  the  conditions  of  interest  in 
most  condensation  problems  the  rotational  mode  of  mole- 
cules may  be  taken  to  be  in  equilibrium  and  vibrational 
modes*  are  unexcited.  Fortunately  therefore,  only  one 
active  relaxation  mode  needs  be  considered  in  homogeneous 
nucleation  problems.  Thus,  the  mass  fraction  of  conden- 
sate may  conveniently  be  taken  as  the  progress  variable, 

<j> , in  Eq.  (3.2).  However,  the  formalism  of  a general 
rate  process  applies,  and  a relaxation  equation  of  the 
general  form: 


a£  = *(p,  s,<i>)  , 


(3.3) 


arises,  where  f is  a prescribed  function  of  p,  s,  and  <J> . 
For  small  deviations  from  equilibrium,  Eq.  (3.3)  may  be 
linearized  to  give 


d*  = _ (*-v 

dt  x 


(3.4) 


Here,  $ is  the  equilibrium  value  of  the  progress  variable 
<})  (e.g.  the  condensate  mass  fraction)  under  the  local  con- 
ditions, and  t = T(p,s,<f>)>0  is  the  relaxation  time  with  the 
property  that  at  equilibrium, 


t ( p , s ,<f>  ) = 0. 
e 


(3.5) 


* This  latter  mode  may  have  to  be  considered  in  conden- 
sation studies  of  polyatomic  molecules. 


Interesting  work  in  gasdynamics  concerns  the  behavior 
of  these  equations  and  it  is  discussed  in  different 


recent  books . * 


More  recent  reviews  by 


Marble  (1969)  on  two-phase  flows,  Becker  (1972)  on 
chemically  reacting  flows,  and  by  Wegener  (1975)  on  condensing 
nonequilibrium  flows  are  also  of  interest  in  this  connection. 

One  additional  fundamental  problem  that  arises  is  the 
proper  definition  of  the  speed  of  sound  in  a condensing 
medium.**  In  a non-reacting  medium,  small  disturbances 
propagate  at  the  local  acoustic  velocity  defined  unambiguous- 
ly by  Eq.  (2.10).  Here,  the  differential  in  Eq.  (2.10) 
is  well  defined  since  p may  be  expressed  as  a function 
of  p and  s using  Eq.  (3.1).  In  a condensing  mixture, 
however,  p is  a function  of  three  (or  more)  variables, 


p = p(p  , s ,4> ) , 


(3.6) 


as  in  Eq.  (3.2),  and  thus  there  is  an  infinite  number  of 
differentials  having  the  form  Op/9p)s.  Indeed  in  con- 
densing flows  we  deal  with  the  acoustics  of  dispersive 
media  in  treating  the  propagation  of  small  disturbances. 


* In  addition  to  books  cited  at  the  beginning  of  Section 
2,  see  Clarke,  J.F.  and  McChesney,  M. , Dynamics  of  Real 
Gases,  Butterworth,  London  (19  64),  and  Wegener,  P.P.  CEd . ) 
Nonequilibrium  Flows  Parts  I and  II,  Marcel  Dekker,  N.Y. 
(1969  , 197  0')“ 

**  Again  the  term  condensing  medium  is  to  be  interpreted 
broadly  to  include  systems  involving  chemical  reaction 
and  other  finite  rate  processes. 
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There  are  two  limits  with  important  physical  signifi- 
cance in  taking  the  derivative  of  the  right  hand  side 
of  Eq . (3.6): 

a,  = A3'p/3p)  7 , (3.7) 

± S 5<P 

and 

a = /( 3p/ 3 p ) " ~ . (3.8) 

e s,  <})  = (t>e 

Here,  a^  and  ag  are  defined  as  the  frozen  and  equilibrium 
sound  speeds,  respectively.  The  differentiation  in  Eq. 
(3.7)  is  taken  with  <p  held  constant  regardless  of  the 
state  of  the  gas.  Thus  the  composition  of  the  mixture 
is  fixed,  or  "frozen".  On  the  other  hand,  the  differenti- 
ation in  Eq.  (3.8)  is  taken  with  <J>  retaining  its  equilib- 
rium value  appropriate  for  the  prevailing  state  of  the  gas. 
The  equilibrium  condition 

4>e  = 4>e  (p,s)  , (3.9) 

is  given  by  the  law  of  mass  action  for  chemical  reactions, 
the  Clausius-Clapeyron  equation  for  condensation,  or 
equivalents  for  other  processes.  When  Eq.  (3.9)  is  in- 
troduced in  Eqs.  (3.6)  and  (3.8),  the  classical  gasdynamic 
result,  Eq.  (2.10)  is  recovered  as  expected.  The  terms 
"frozen"  and  "equilibrium"  may,  in  turn,  be  interpreted 
to  indicate  the  limiting  behavior  of  the  flow  when  the 
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rates  of  condensation  are  infinitely  slow  and  infinitely 
rapid,  respectively.  We  then  have  t = 00  or  t = 0 in  Eq. 
(3.4)  . 

Chu  (1958)  and  Broer  (1951)  both  showed  that 
the  characteristic  speed  of  propagation  of  small  dis- 
turbances in  a reacting  medium  is  the  local  frozen  sound 
speed  in  all  reacting  flows  having  a finite  reaction 
rate.  Only  when  the  reaction  rate  is  truly  infinite, 
i.e.,  when  the  flow  is  exactly  in  thermodynamics  equi- 
librium, is  the  characteristic  speed  given  by  the  equi- 
librium sound  speed.  Thus,  the  variation  of  character- 
istic speed  with  reaction  rate  is  discontinuous,  being 
a^  at  all  finite  rates  and  ag  only  when  the  rate  is  in- 
finite (Chu  1958). 
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The  definition  of  characteristic  speed  in  a react- 
ing flow  is  important  because  it  serves  to  define  the 
proper  sound  speed,  the  Mach  number,  and  the  direction 
of  characteristics  for  the  partial  differential  equa- 
tions describing  the  flow.  It  can  be  shown  using  thermo- 
dynamic stability  arguments  that 

a^  > a , (3.10) 

f e 

and  typically 

1.0  < a .p/a  < 1.2  (3.11) 

f e 

for  dissociating  gases.  For  condensing  flows  we  find: 

1.1  < ap/a  < 1.8  (3.12) 

f e 

for  H2O  (Wegener  and  Mack  1958)  resulting  in  a much  more 
pronounced  effect  of  condensation  on  the  dispersion  of 
sound  than  the  effect  found  in  reactions. 

Nozzle  flow  of  gases  undergoing  dissociation  or  other 
single  chemical  reactions  is  well  understood  at  this  time. 
When  the  flow  starts  at  a high  temperature  in  the  supply 
in  thermodynamic  equilibrium  the  gas  may  be  partly  disso- 
ciated. Recombination  sets  in  as  the  gas  mixture  is  cooled 
in  the  expansion.  At  some  point  determined  largely  by  the 
local  collision  frequency  of  the  gas,  "sudden  freezing" 
occurs  as  predicted  by  Bray  (1959),  the  reaction  ceases 
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and  the  flow  departs  from  equilibrium.  An  interesting 
parallel  may  be  drawn  for  nozzle  flows  with  condensa- 
tion (Wegener  1975).*  In  vapor  expansions  condensation 
is  initially  absent  in  the  nozzle  although  the  vapor 
may  be  highly  supersaturated  and  the  mixture  composition 
remains  constant.  Hence  in  contrast  to  dissociating 
flow  the  expansion  starts  in  the  "frozen"  condition. 

With  the  onset  of  condensation  by  nucleation  and  droplet 
growth  saturation  is  attained  quickly  resulting  in  rapid 
vapor  depletion.  Henceforth  the  remaining  vapor  in  the 
gas  phase  in  the  flow  stays  in  equilibrium  with  the 
finely  dispersed  liquid  phase.  Thus  condensation  occurs 
spontaneously  by  "sudden  equilibrium"  in  an  exactly 
opposite  process  to  that  of  "sudden  freezing"  discussed 
above  for  reacting  flow.  Analytical  work  has  recently 
been  done  by  Petty  et  al.  (1972)  in  an  attempt  to  locate 
the  "sudden  equilibrium"  onset  point  of  condensation 
through  some  simple  criteria  in  the  spirit  of  Bray's  work. 

In  sum,  we  note  that  much  experimental  and  theoretical 
work  has  been  done  on  nonequilibrium  nozzle  flows  and  recent 
advances  in  the  field  have 

* This  suggestion  was  made  by  Prof.  N.H.  Johannesen  of 
Manchester  University  in  a discussion  of  an  earlier  paper 
by  Wegener  and  Parlange  (1967)  on  p.  633. 
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been  reviewed  by  Bray  (1970).  Steady,  shockless  flows 
in  slender  nozzles  involving  a single  nonequilibrium 
process  are  well  understood  at  this  time.  This 

knowledge  has  been  successfully  used  in  recent  years  to 
determine  the  rate  constants  of  some  chemical  processes 
through  gasdynamic  measurements  of  the  nozzle  flow  and 
also  to  determine  the  rate  of  condensation  by  homogeneous 
and  binary  nucleation  and  cluster  growth, 
b)  Flow  with  Heat  Addition 

The  energy  equation  for  one-dimensional , adiabatic 
and  frictionless  flow,  Eq.  (2.19)  may  easily  be  modified 
to  admit  heat  exchange  with  the  surroundings  resulting 
in  "diabatic  flow."  We  shall  not  presently  be  concerned 
with  the  kinetics  of  this  heat  transfer;  rather,  we  shall 
focus  on  the  gasdynamic  effects  on  the  flow  once  the  heat 
transfer  has  been  completed.  This  view  is  consistent 
with  our  previous  assumption  of  high  Reynolds  number  flow 
with  thin  boundary  layers.  In  compressible  flow  of  real 
fluids,  there  will,  in  general,  be  temperature  gradients 
along  and  normal  to  the  streamlines,  as  dictated  by  the 
conditions  at  the  flow  boundary  (e.g.,  heat  conducting  or 
insulating  sidewalls  held  at  fixed  or  variable  temperature). 
Like  velocity,  it  is  found  experimentally  that  temperature 
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variations  across  the  streamlines  are  also  restricted 
to  a narrow  region  near  the  wall.  In  fact,  a thermal 
boundary  layer  exists  near  the  wall.  Outside  this 
layer  the  heat  conductivity  of  the  fluid  plays  no  role, 
and  the  fluid  behaves  as  though  it  were  nonconducting. 
The  ratio  of  the  thickness  of  the  thermal  and  the 
viscous  boundary  layers  respectively  is  proportional 
to  the  inverse  square  root  of  the  so-called  Prandtl 
number.  This  dimensionless  parameter  is  defined  by: 


Pr 


nc^ 

K 


(3.13) 


with  k,  the  thermal  conductivity.  The  Prandtl  number 
has  the  physical  significance  of  expressing  the  ratio 
of  the  diffusivity  of  viscous  to  thermal  effects  res- 
pectively and  for  most  gases  we  find  Pr^l . Therefore, 
the  thermal  and  viscous  boundary  layers  are  roughly  of 
the  same  thickness.  In  fact,  in  the  one-dimensional 
treatment  of  nozzle  flows  using  the  effective  area 
distribution  given  by  a static  pressure  calibration 
discussed  previously,  both  the  viscous  and  the  thermal 
boundary  layers  are  accounted  for,  and  the  flow  may  be 
treated  as  inviscid  and  nonconducting.  Within  the 
boundary  layer,  the  temperature  increases  from  that 
given  by  the  isentropic  relation  Eq.  (2.22)  roughly  to 
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the  stagnation  temperature,  Tq  , at  the  wall  for  those 
gases  with  Pr«l. 

In  the  following  discussion,  therefore,  we  can 
assume  that  the  gas  is  nonconducting,  and  heat  addi- 
tion in  one-dimensional  flow  is  uniform  at  a given 
cross  section  normal  to  the  flow  direction.  The 
amount  of  heat  added  to  unit  mass  of  the  flow,  q may, 
of  course,  vary  with  x and  the  function  q(x)  depends 
on  the  process  at  hand.  Heat  addition  of  this  kind 
is  realized  in  flows  subject  to  thermal  radiation, 
flows  with  chemical  reaction*  (e.g.  combustion),  flows 
with  condensation*,  etc. 

In  the  diabatic  flows  discussed,  the  continuity, 
momentum  and  thermodynamic  state  equations,  (2.14), 
(2.15),  and  (2.16),  remain  valid.  The  energy  equation, 
Eq.  (2.19),  for  steady,  adiabatic  flow  of  an  inviscid 
gas  must  now  be  modified  by  including  a term  dq  to  admit 
external  heat  exchange  , 

dh  + udu  = dq.  (3.14) 


* It  may  be  argued  that  in  these  cases  no  heat  addition 
from  an  external  source  takes  place.  Heat  is  merely 
being  transformed  from  one  form  to  another.  However,  it 
is  often  more  convenient  and  informative  to  consider 
flows  with  chemical  reactions  or  condensation  external 
heat-addition  problems  by  a suitable  choice  of  the  system. 
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For  an  ideal  gas , it  is  convenient  to  regard  the  heat 
addition  as  a change  in  the  local  stagnation  tempera- 
ture , * i . e.  , 

c_dT  + udu  = cDdT  , (3.15) 

P P o 

with 

dq  = cndT  . (3.16) 

^ P o 


Detailed  discussions  of  the  influences  of  various  pro- 
cesses on  the  gasdynamics  variables,  M,  u,  p,  p,  etc., 
of  one-dimensional  flow  are  available  in  Shapiro's 
textbook.**  For  example,  in  diabatic  flow  the  changes 
in  Mach  number,  cross-sectional  area,  and  stagnation 
temperature  are  related  by 


dMa  _ 2(1+1^M2)  dA  (i+YM2)(l+^iM2)  dTQ 


M 


1-M‘ 


1-M2 


(3.17) 


It  can  be  seen  from  Eq . (3.17)  that  heat  addition  (dTQ>0) 
to  subsonic  flow  (M<1)  is  accompanied  by  an  increase  in 


* With  heat  addition,  T ^ T where  T is  e.g.  the 

o o i o 1 

temperature  in  the  supply  of  a nozzle.  Often  the  sub- 
script 1 is  dropped,  however. 

**  Shapiro,  A.H.,  The  Dynamics  and  Thermodynamics  of 
Compressible  Fluid  Flow,  Vols.  I and  II,  Ronald  Press, 
N.Y.  (1953). 


Mach  number,  while  that  to  supersonic  flow  (M>1)  de- 
creases the  Mach  number.  Thus,  heat  addition  always 
tends  to  drive  the  flow  Mach  number  towards  unity. 
Moreover,  since  the  coefficients  for  dA  and  dT  in  Eq . 
(3.17)  are  of  opposite  sign,  the  effects  of  heat  addition 
on  M are  equivalent  to  and  may  be  duplicated  by  an  area 
reduction.  This  qualitative  analogy  between  heat  addi- 
tion and  area  reduction  is  found  to  be  generally  valid 
for  the  influence  on  the  other  flow  variables  as  well. 

Using  the  area  reduction  analogy  to  heat  addition 
we  note  that  for  a given  Mach  number  there  is  a specific 
amount  of  heat  that  may  be  added  to  the  gas  stream  which 
is  equivalent  to  a reduction  of  the  cross-sectional  area 
to  the  critical  throat  area  A*.  The  Mach  number  after 
this  critical  heat  addition  will  then  be  unity  and  at 
this  point  the  flow  becomes  "thermally  choked".  This 
maximum  or  critical  amount  of  heat  addition,  qm=  . (for 
condensation  note  Lukasiewicz  and  Royle  1948),  is  given 
by 

qmax  _ (M2  -1 ) 2 . (3.18) 

cpTo  M2(y2-1)[M2+2/(Y-1)] 

Heat  addition  exceeding  the  value  of  q may  cause 
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shock  waves  or  unsteady  flows  to  occur  (e.g.  Wegener 
1975),  and  these  phenomena  invalidate  the  simple  treat- 
ment of  Section  2.  In  this  work, we  shall  only  be  in- 
terested in  diabatic  flow  where  the  heat  addition  does 
not  lead  to  thermal  choking  or  discontinuities  in  the 
flow.  Therefore  , the  basic  conservation  equations  in 
the  form  of  Eqs . (2.14),  (2.15)  and  (3.14)  remain  valid. 

: 

I 

3c.  Flow  with  Condensation 

Condensing  flows  show  simultaneous  nonequilibrium 
and  heat  addition  effects  and  the  progress  variable  <j> 
for  Eqs.  (3.2)  and  (3.6)  may  be  identified  with  the 
mass  fraction  g of  the  condensate.  The  thermodynamic 
state  of  the  flow  is  fixed  once  g and  two  other  thermo- 
dynamic variables,  say  p and  T,  are  specified.  The 
incremental  heat  addition  per  unit  mass  of  mixture  of 
the  flow  is  related  to  the  condensate  mass  fraction  by: 

dq  = L^dg.  (3.19) 

Here  ,L  stands  for  the  latent  heat  of  the  bulk  conden- 
* 00 


sate  and  possible  curvature  effects  on  this  value 
(Wegener  and  Parlange  1967)  are  not  considered. 
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Provided  that  the  pressure  ratio  across  the  nozzle 
is  sufficient  to  produce  sonic  flow  at  the  nozzle 
throat,  and  that  supercritical  heat  addition  does  not 
take  place,  integration  of  the  gasdynamic  equations  for 
supersonic  nozzle  flows  are  rather  straightforward. 

The  theory  of  supersonic  nozzle  flow  with  vapor  conden- 
sation was  first  developed  by  Oswatitsch  (1942).  He 
showed  that  the  presence,  in  the  flow,  of  condensate 
droplets  produced  by  homogeneous  nucleation  does  not 
invalidate  the  assumption  of  inviscid  flow  because 
they  are  usually  extremely  small  (typically  r £ 1000,8) 
and  thus  they  move  practically  with  the  gas  stream.* 
Hence,  the  viscous  drag  between  the  droplets  and  the 
gas  stream  may  be  neglected  outside  the  region  of  extreme 


* For  proof  of  this  statement  see  e.g.  Wegener  (1969). 
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velocity  gradients  as  found  in  shock  waves  or  boundary 
layers.  It  may  be  argued  that  although  the  momentum 
transport  between  the  gas  stream  and  the  droplets  is 
negligible,  the  heat  transfer  across  the  phase 
boundary  may  not  be  ignored,  and  heat  conduction  ought  to 
be  included  in  the  energy  equation.  In  fact,  it  is  due 
to  the  transfer  of  latent  heat  of  condensation  from  the 
droplets  to  the  gas  phase  that  condensation  is  "felt" 
gasdynamically . However,  since  the  condensate  is  found 
to  be  finely  dispersed  (~1012/cm3)  and  therefore  the  droplets 
may  be  considered  as  uniformly  distributed  heat  sources 
throughout  the  gas  stream.  Thus,  the  heat  may  be  viewed 
as  being  generated  internally  in  the  gas  phase  and  the 
problem  of  heat  transfer  between  phases  may  be  circum- 
vented for  the  time  being.* 

Oswatitsch’s  assumption  that  condensation  proceeds  by 

homogeneous  nucleation  and  droplet  growth  enabled  him  to 

compute  the  mass  fraction  of  condensate  as  a function  of 

the  position  along  the  nozzle  by  summing  up  the  mass  of  all 
droplets 

* On  a different  length  scale,  for  example,  in  distances 
on  the  order  of  the  diameter  of  the  droplet,  this  heat 
transfer  process  may  be  important.  In  fact  it  is  one 
limiting  factor  of  the  droplet  growth  rate. 
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present.  He  realized  that  although  droplets  of  dif- 
ferent sizes  would  be  present  at  any  given  location 
along  the  nozzle,  it  was  possible  to 
determine  the  size  of  each  individual  droplet  by 
tracing  its  history  back  to  its  origin  as  a nucleus 
of  critical  size  and  by  integrating  the  growth  rate 
along  the  particle  path.  Furthermore,  he  assumed  that 
no  disintegration  or  coalescence  of  droplets  occurs 
resulting  in  the  conservation  of  the  number  of  droplets. 
Thus,  both  the  total  number  of  droplets  and  the  size 
of  each  may  be  determined  yielding  the  mass  fraction 
of  condensate.  A condensation  rate  equation,  in  the 
general  form  of  Eq.  (3.3)  may  therefore  be  derived  from 
a known  droplet  growth  rate  and  nucleation  rate.  For 
this  purpose,  Oswatitsch  (1942)  gave  a droplet  growth 
equation,  and  he  used  the  steady  state  homogeneous 
nucleation  rate  equation  of  Becker  and  DcJring 
(1935)  to  describe  the  kinetics  of  condensation.  The 
droplet  growth  law  was  based  on  an  approximate  free- 
molecule  treatment  of  the  simultaneous  heat  and  mass 
transfer  occurring  at  the  surface  of  the  droplet. 

For  this  assumption  to  be  valid,  the  droplet  diameter 
must  be  much  smaller  than  the  mean  free  path  of 
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the  surrounding  gas,  a condition  satisfied  in  flows  of 
vapor  in  the  presence  of  an  inert  carrier  gas.  Typi- 
cally for  moist  air,  where  the  maximum  droplet  radii 

are  about  100  X at  the  completion  of  condensation  and 

. o 

the  mean  free  path  is  of  the  order  of  1000  A,  Eq.  (2.3) 

shows  Kn~10.  For  condensation  of  a pure  vapor  (e.g., 

steam) , the  droplets  in  the  condensation  zone  eventu- 

o 

ally  grow  beyond  1000  A.  However,  in  the  region  immedi- 
ately downstream  of  the  location  where  the  droplets  are 
formed,  they  are  again  of  a sufficiently  small  size  for 
the  free  molecule  flow  assumption  to  hold. 

Under  these  assumptions,  the  primary  effect  of  con- 
densation on  the  supersonic  nozzle  flow  is  the  addition 
of  the  latent  heat  of  condensation  to  the  gas  stream.  Therefore, 
the  gasdynamic  equations  of  motion  are  similar  to  those  of 
one-dimensional,  frictionless,  diabatic,  and  steady  flow, 
discussed  before.  If  condensation  occurs  only  in  the 
supersonic  section  of  the  nozzle,  integration  in  the  flow 
direction  may  conveniently  start  at  the  throat,  using  the 
sonic  conditions  as  a starting  point.  This  situation 
avoids  integration  through  the  nozzle  throat,  which  is  a 
mathematical  singularity  as  noted  in  the  cited  textbooks. 

Thus,  Oswatitsch  was  able  to  perform  a stepwise  numer- 
ical integration  of  the  system  of  equations  for  steady 


supersonic  nozzle  flow.  In  the  pre-computer  days  of 
1942  he  obtained  remarkable  agreement  of  the  theory 
with  steam  experiments  (Yellot  and  Holland  (1937), 
Binnie  and  Woods  (1938)), and  his  own  experiments  with 
moist  air.  In  the  intervening  years,  Oswatitsch's 
formulation  was  consistently  proven  to  be  successful  as 
adapted  by  many  later  investigators.  Moreover, 
Oswatitsch’s  model  has  conversely  become  the  basis  for 
using  supersonic  nozzles  in  the  experimental  study  of 
condensation  phenomena.* 

Along  similar  lines,  one  of  us  (B.J.C.  Wu)  has 
derived  equations  for  the  one-dimensional  unsteady 
expansions 


* Papers  using  this 

approach  include  Wegener  and  Mack  (1958),  Stever  (1958), 
Pouring  (1963),  Hill,  Witting  and  Demetri  (1963), 
Gyarmathy  and  Meyer  (1965),  Griffin  and  Sherman  (1965), 
Duff  and  Hill  (1966),  Hill  (1966),  Stein  (1967),  Wegener 
and  Parlange  (1967),  Deych  and  Filipov  (1968),  Dawson 
et  al.  (1969),  Jaeger  et  al.  (1969)  Puzyrewski  (1969), 
Wegener  (1969),  Roberts  (1969),  Chmielewski  and  Sherman 
(1970),  Clumpner  (1970),  Barschdorff  (1970,  1971),  Wu 
(1972),  Wegener,  Clumpner  and  Wu  (1972),  Barschdorff  et 
al.  (1972),  and  Wu  (1974). 
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with  condensation  associated  with  shock  tube  flows  valid 
when  the  latent  heat  added  to  the  gas  stream  is  rela- 
tively small  as  used  by  Barschdorff  (1975)  and  as  seen 
in  Section  6 . 

4 FLOW  SYSTEMS  FOR  CONDENSATION  RESEARCH, 
a)  Steady  Nozzle  Flow  and  Ludwieg  Tube. 

Nozzles  provided  the  first  data  on  steam  condensa- 
tion (e. g. ,Stodola  1927):  they  were  used  by  Oswatitsch 
(1942)  for  his  classical  studies  of  nucleation,  and 
have  since  been  applied  in  many  laboratories.  Such 
Laval*  nozzles  are  most  useful  in  condensation  research 
if  the  contour  is  designed  for  continuously  expanding 
flow  (e.g.  Wegener  1964)  rather  than  for  a uniform  Mach 
number  distribution  in  the  test  section  as  found  in 
supersonic  wind  tunnels.  In  the  expansion,  condensa- 
tion causes  a deviation  of  the  pressure  distribution 
from  that  observed  in  isentropic  flow  outside  the  bound- 
ary layers  as  determined  from  pressure  surveys  (see 
Section  2).  Examples  of  such  static  pressure  traces, 
p(x)  , as  obtained  by  a continuous  pressure-distance 


* The  Swedish  engineer  de  Laval  introduced  converging- 
diverging  nozzles  (Section  2b)  in  steam  turbines  in  1883. 


survey  in  two  different  nozzles  are  shown  in  Fig.  4.1. 

The  duration  of  the  steady-state  in  such  experiments  is 
of  lesser  importance  and  continuous  as  well  as  intermit- 
tent wind  tunnels  are  used.  The  choice  of  the  mode  of 
operation  depends  on  the  capacity  of  the  vacuum  equipment, 
high  pressure  supply,  etc.,  and  on  the  cost  of  the  gas 
mixture  to  be  expanded.  In  any  case,  the  maximum  expan- 
sion ratio  of  the  nozzle,  A/A*,  or  the  corresponding 
Mach  number  must  be  chosen  to  be  beyond  the  initially 
estimated  point  of  departure  from  the  isentropic  expansion  as 
shown  in  Figs.  4.1  and  4.2.  The  nozzle  is  then  designed 
with  the  aid  of  Eq.  (2.31)  and  the  corresponding  values 
of  the  flow  table  in  Appendix  1. 

An  advantageous  application  of  a nozzle  to  conden- 
sation studies  is  provided  by  the  so-called  Ludwieg  tube 
(Ludwieg  1955),  a combination  shock  tube  and  nozzle 
facility.  This  device  provides  steady  flows  in  the  milli- 
second to  second  range  (depending  on  the  length  of  the  supply 
tube)  while  operating  without  valves  and  other  ma- 

chinery. Here, non-equilibrium  flows  of  poisonous  or 
otherwise  dangerous,  or  costly,  gases  can  be  produced 
since  only  a small  anount  of  a carefully  prepared  mixture 
is  required.  The  instrumentation  for  condensation 
studies  (Wegener  and  Cagliostro  1973) 
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requires  pulse  techniques,  however,  it  is  in  the 
range  of  common  practice  in  experimental  gasdynamics  . 

The  application  of  the  gasdynamic  equations  for 
diabatic  flow  (Section  3b)  to  experiments  such  as  those 
of  Fig.  4.1  leads  to  the  results  of  Fig.  4.2.  Ordinarily 
the  assumption  of  dealing  with  a thermally  and  calori- 
cally  perfect  gas  can  be  made  and  we  obtain  all  local 
state  parameters  from  a solution  of  the  equations  of 
motion  in  conjunction  with  known  A(x)  and  measured  p(x). 

In  our  example  the  heat  addition  due  to  condensation  started 
at  x = 1.6  cm  and  it  is  shown  by  the  computed  increase  of 
the  local  total  (stagnation)  temperature  above  its  value 
in  the  supply  whose  properties  carry  the  subscript  oi , 
to  distinguish  it  from  the  computed  local  stagnation  pro- 
perties. Correspondingly,  the  irreversible  character  of 
the  heat  addition  process  is  noted  in  a drop  of  stagnation 
pressure.  It  is  important  to  recall  that  up  to  this  point 
no  assumptions  on  the  condensation  process  were  required 
and  the  thermodynamic  state  throughout  the  condensation 
zone  could  be  computed  from  a simple  experiment.  Only 
after  later  applying  Eq . (3.19)  in  Section  6 to  data  such 
as  those  of  Fig.  4.2  shall  we  arrive  at  a determination 
of  the  condensate  mass  fraction. 
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Among  the  additional  diagnostic  methods, 

light  scattering  offers  an  attrac- 
tive tool.  In  Fig.  4.3  from  Stein  (1969a)  we  see  results 
for  a cloud  of  water  droplets  obtained  by  laser  light 
scattering  beyond  the  nozzle  exit  of  the  equipment  des- 
cribed in  the  reference  and  elsewhere  (Stein  1969b). 
Similar  data  were  obtained  by  Roberts  (1969)  and  Clumpner 
(1971)  who  worked  with  ethanol.  Since  the  wavelength 
of  the  laser  X>>r,  and  multiple  scattering  is  negligible, 
we  are  dealing  with  pure  Rayleigh  scattering,  for  which 
we  have 


I 

o 


h&2  16tt  4 
R2  A4 


cos20 


m2-I  2 
m2  +2 


f(r)dr  , 


(4.1) 


where  I and  IQ  are  the  intensities  of  the  incident  and 
o u 

scattered  radiation  respectively  in  units  of  energy /time. 
Furthermore,  h and  Jt  are  the  dimensions  of  the  viewing 
window  at  angle  0 at  a distance  R from  the  scattering 
volume.  The  index  of  refraction  of  the  clusters  is  given 
by  m,  N is  the  total  number  of  clusters  per  unit  volume, 
and  f(r)  is  the  normalized  distribution  function.  Eq. 
(4.1),  is  valid  for  the  case  where  the  electric  vector 
of  the  plane-polarized  laser  light  is  parallel  to  the 
plane  of  observation  (E^  ) , while  for  the  case  with  the 
electric  vector  perpendicular  (E^) , the  cos20  factor 
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would  not  appear.  Fortunately , a second  equation  is  at 
hand  to  link  N and  r as  required  to  solve  Eq . (4.1). 
Since  we  know  (e.g.  Wegener  andiPouring  1964)  that 
practically  all  the  mass  of  the  condensing  vapor  con- 
denses in  the  condensation  zone  upon  reaching  the  equi- 
librium state*  at  the  observation  station  used  in  Fig. 
4.3,  we  can  write  for  the  conservation  of  mass 


oo 

U)qPt  3 ^ ^ r3  f (r )dr  . 


(4.2) 


Different  distribution  functions  for  the  radii  can  be 
matched  to  the  light  scattering  signal  in  Eqs.  (4.1)  and 
(4.2).  Figure  4.3,  however,  gives  a simple  evaluation 
assuming  a delta  function  size  distribution,  i.e.  a 
uniform  radius  for  all  clusters  (see  Section  6). 

Finally,  density  distributions  by  interferometry, 
probing  wi.th  pitot  tubes,  etc.  have  all  been  used  in 
nozzles  and  we  refer  to  such  techniques  in  the  list  of 
papers  in  Appendix  2 . 

All  these  techniques  enable  us  to  determine  the 
pressure  and  temperature  on  the  isentrope  at  which  the 
initial  condensate  appearance  is  noted. 


* For  proof  see  Section  6. 
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This  "onset"  of  condensation  yields  the  adiabatic  super- 
cooling, AT^,  of  Section  1,  and  Fig.  1.2,  or  with  known 
p^CT),  it  permits  us  to  compare  the  critical  supersatura 
tion  (p/p^)*  customarily  employed  in  cloud  chamber 
research.  Figure  4.4  gives  such  "raw  data"  on  water 
vapor  condensation  in  steady  nozzle  flow  for  steam  or 
in  carrier  gases.  The  sources  of  these  data  are  listed 
in  Table  4.1.  The  scatter  of  the  points  can  be  ratio- 
nalized by  remembering  Eq.  (1.1)  and  Fig.  1.4  since 
different  cooling  rates  and  different  initial  conditions 
are  shown.  However,  a "Wilson  line"  shifted  with  re- 
spect to  the  equilibrium  vapor  pressure,  p^CT),  emerges. 

b)  Shock  Tube. 

There  are  several  useful  applications  of  the  shock 
tube  involving  different  regions  of  the  flow  (Fig.  2.2). 
Historically  the  first  application  of  the  shock  tube  to 
condensation  studies  (Wegener  and  Lundquist  1951)  re- 
lates to  the  expansion  fan  between  regions  4 and  3 in 
which  the  gas-vapor  mixture  undergoes  an  unsteady,  isen- 
tropic  expansion  governed  by  the  equations  of  Section  2c 
There  is  some  question  on  the  quality  of  the  flow  set 
in  motion  by  the  fan.  However, 
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experimental  pCt)  histories  agree  perfectly  with  the 
prediction  of  the  ideal  theory  provided  an  initial 
shift  of  the  spatial  coordinate  has  been  made  to  account 
for  inertial  and  other  effects  related  to  diaphragm  rupture.* 

Flow  fields  were  photographed  by  Glass  and  Patterson  (1955) 
who  show  beautiful  experimental  x-t  diagrams  with  condensation 
fronts  that  are  the  unsteady  analog  to  our  steady 
nozzle  flow  photograph  Fig.  1.1.  Pressure-time  his- 
tories without  and  with  condensation, respectively,  appear 
like  the  p(x)  graphs  in  Fig.  4.1  (Wegener  1975, 

Barschdorff  1975).  At  a given  observation  station,  a pressure 
rise  with  respect  to  the  ideal  isentropic  expansion  and/or 
a light  scattering  signal  reveal  the  onset  of  condensa- 
tion and  the  history  of  the  corresponding  particle  path 
(Fig.  2.2)  can  be  computed  as  shown  in  Section  2. 

Usually  the  cooling  rate  is  lower  in  such  fans  than  in 
a nozzle.  In  fact,  at  a sufficient  upstream  distance 
from  the  diaphragm,  the  low  cooling  rates  of  cloud 
chambers  may  be  duplicated.  Finally,  in  order  to  avoid 
any  effects  of  heat  addition  from  condensation  that  may 
spread  in  both  directions  in  the  subsonic  outflow  of  the 

* This  empirical  correction  can  be  avoided  by  simult- 
taneous  observations  at  two  or  more  stations  (Lee  1974). 
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driver  section*,  the  onset  of  condensation  may  be 
arranged  to  coincide  with  the  tail  of  the  expansion 
fan  as  shown  in  Fig.  2.2  (Barschdorff  1975).  Since 
it  is  normally  impractical  to  shift  the  observation 
station,  this  condition  can  be  achieved  by  a judi- 
cious choice  of  p^/p^  ^or  Eq.  (2.44).  Onset  data 
obtained  in  shock  tubes  may  then  be  treated  just  like 
those  for  nozzles. 

Other  schemes  are  useful  for  condensation  studies 
of  materials  such  as  metals  that  are  vapors  only  at 
high  temperatures.  Placed  in  the  driven  section  the 
substance  of  interest  will  first  be  heated  and  va- 
porized by  the  shock  wave.  Kung  and  Bauer  (1971) 
have  thus  generated  iron  vapor  at  T~  2,600  K.  Making 
use  of  earlier  ideas  on  the  so-called  tailored  inter- 
face technique  (e.g.  Wittliff  et  al . 1959),  they 
produced  controlled  expansions  in  which  to  study  the 
condensation  by  homogeneous  nucleation  of  the  previously 


* Note  that  the  tail  of  the  expansion  fan  moves  to  the 
left  in  Fig.  2.2,  i.e.  u^ca^. 


vaporized  iron  vapor.  However,  these  results  as  well 
as  the  use  of  the  shock  heated  region  2 itself  for 
the  precipitation  of  e.g.  lead  held  at  constant  tem- 
perature (Homer  et  al.  1971,  Homer  and  Hurle  1972), 
have  not  yet  led  to  satisfactory  comparisons  with 
nucleation  theory. 

c)  Molecular  Beam. 

Becker  and  Henkes  (1956)  discovered  that  rapid 
expansions  in  high  intensity  molecular  beams  often 
lead  to  clustering  of  various  species.  Such  beams 
are  produced  by  expanding  gases  through  hole-type 
nozzles  into  vacuum  leading  to  very  high  Mach  number 
flows  in  the  free  molecular  regime  (Kn>>l).  The 
cooling  rates  in  such  expansions  range  up  to  109  C/sec, 
and  it  has  been  shown  that  the  steady-state  assumption 
of  classical  nucleation  theory  is  no  longer  valid 
under  these  circumstances  (Andres  and  Boudart  1965). 

In  spite  of  the  fasinating  cluster  studies  that  this 
technique  permits,  we  shall  omit  a discussion  in  this 
chapter.  For  one,  our  primary  aim  to  compare  well  con- 
trolled gasdynamic  continuum  data  with  classical  theory 
breaks  down  for  the  above  reasons.  Furthermore, 

Hagena  (1974)  has  just  published  an  excellent  review 


of  the  current  status  of  condensation  in  molecular 
beams  to  which  we  shall  defer. 

5 EQUATIONS  OF  MOTION  WITH  HOMOGENEOUS  NUCLEATION 
AND  DROPLET  GROWTH. 

a ) Gas  dy nami cs . 

In  this  section,  we  shall  derive  the  governing 
equations  for  supersonic  nozzle  flow  with  condensa- 
tion.* These  equations  will  be  derived  under  the 
assumptions  laid  out  in  Sections  2 and  3,  and  we  shall 
focus  on  the  special  features  of  flows  with  conden- 
sation . 

For  the  gaseous  mixture  (excluding  the  droplets 
of  condensate),  the  equation  of  state,  Eq.  (2.5),  for  a 
perfect  gas  is  assumed.  Here,  the  unscripted  symbols 
refer  to  the  properties  of  the  gaseous  mixture,  while 
the  properties  of  the  two-phase  mixture  are  denoted 
by  the  subscript  T.  Since  the  condensate  occupies  a 
negligible  volume  compared  with  that  of  the  gas  phase, 
we  may  write 

pT  = p + pc  = Pf  + pv  + pc  ’ (5,1) 

We  recall  that  is  the  condensate  density  referred  to 
the  volume  of  the  gas  phase,  and  is  the  (apparent) 

* Similar  derivations,  using  slightly  different  notations, 
may  be  found  in  Wegener  and  Mack  (1958),  Stever  (1958), 
Hill  (1966),  and  Wegener  (1969). 
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total  density  of  the  two-phase  medium.  This  total  den- 
sity is,  of  course,  meaningful  only  on  length  scales 
large  compared  with  the  average  inter-droplet  distance. 
With  g,  the  mass  fraction  of  condensate,  we  have: 


and 


PA  = gPn 


pT  = p/(l-g) 


(5.2) 


(5.3) 


Instead  of  Eq.  (2.8),  the  molecular  weight  of  the  gaseous 
mixture  is  to  be  computed  from : 


1 —co 


o 1 A “o-S  1 


— + 


(5.4) 


p 1-g  pi  1-g  pv 

The  continuity  and  momentum  equations  may  readily  be 
written  in  terms  of  pT , or,  using  Eq.  (5.3): 

puA/ ( 1-g ) = constant  = Q , (5.5) 

and 


_dP  = TT7T  udu  > 


1-g 


(5.6) 


respectively . 


Equation  (5.6)  is  equivalent  to  : 
- Adp  = Qdu  . 


(5.7) 
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The  energy  equation  for  unit  volume  of  the  gas- 
condensate  mixture  is: 

d(pTu2/2  + hT)  = 0 . (5.8) 

A 

Here,  h^,  is  the  enthalpy  per  unit  volume  of  the  two- 
phase  medium.  It  may  be  written  as; 


lm  = p . h . + p h + p'h  , 

T i v v Kc  c ’ 


(5.9) 


where  hu  , hv , and  hc  are  the  enthalpy  per  unit  mass  of 
the  inert  gas,  vapor  and  condensate  respectively.  The 
enthalpies  of  the  vapor  and  condensate, are  related  by 
their  temperature  difference  and  the  latent  heat  of 
condensation.*  Thus: 


h -h  = L + c (T-T  ) , 

VC  C C * 


(5.10) 


where  cc  and  Tc  are  the  specific  heat  and  temperature, 
respectively,  of  the  condensate.  In  most  cases,  we  may 
safely  ignore  the  temperature  difference  between  the  gas 
phase  and  the  condensate  without  introducing  a significant 
error  in  the  energy  equation.  Under  this  assumption, 
we  have 

hy-hc  = L , (5.11) 


* Here  we  shall  again  set  LsL^;  note  Eq.  (3.19). 
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and  the  energy  equation  becomes 

d[pTu2/2  + pihi  + pvhy  + p^(hv-L)]  = 0 . (5.12) 

With  the  definition  of  u>o , the  mass  fraction  of  vapor  in 
the  supply  which  is  assumed  to  be  free  of  condensate,  we 
get 

d[u2/2  + (1-u)  )h.  + id  h -gL]  = 0 . (5.13) 

Defining  cpQ  as  the  constant-pressure  specific  heat  of 
the  vapor-carrier  gas  mixture  in  the  supply,  we  arrive 
at  the  energy  equation  for  flows  with  condensation: 

d[u2 / 2 + cpQT  - gL]  = 0 . (5.14) 

Had  we  not  assumed  equal  temperature  for  the  gas  phase 
and  condensate,  we  would  have  arrived  at  a similar 
equation  with  L replaced  by  L where 

L = L + c (T-T  ) . (5.15) 

c c 

For  a conservative  estimate,  we  figure*  the  error  in  the 
quantity  in  the  square  brackets  in  Eq.  (5.14)  due  to  the 
replacement  of  L by  L to  be  about  0.1%,  which  is  negli- 
gible. Hence,  we  conclude  that  in  the  energy  equation 
for  gas  flow,  we  may  ignore  the  temperature  difference 

* For  typical  conditions  prevailing  in  steam  nozzles, 

T = 300  K,  Tc  = 320  K,  g = 10"2  ,u  = 500  m/sec. 


: 


between  the  droplets  and  the  gas.  The  question  of  drop- 
let temperature  arises  again  in  the  calculation  of 
growth  rates  where  it  will  be  necessary  to  evaluate 
the  droplet  temperature. 


The  presence  of  the  term  gL  in  the  energy  equation,  Eq . 
(5.14)  in  contrast  to  Eq . (2.19)  suggests 
that  vapor  condensation  is  equivalent  to  the  addition 
of  heat  q = gL  to  unit  mass  of  the  flow.  The  equation 
of  state,  Eq.  (2. 6), the  continuity  equation,  Eq.  (5.5), 
the  momentum  equation,  Eq.  (5.6)  or  (5. 7), and  the  energy 
equation,  Eq.  (5.14)  now  complete  the  gasdynamic  des- 
cription of  the  flow  with  condensation.  These  equations 
reduce  to  those  for  perfect  gases  derived  in  Section  2 
when  g is  identically  zero.  Here,  g is  left 

unspecified,  and  we  must  turn  to  the  kinetics  of  the 

condensation  process  to  determine  its  value, 
b)  Condensation  Kinetics. 

The  mass  fraction  of  condensate  g(x)  may  be  computed 

from 


g(x) 


4 IT  ^c 

3 Q 


i 


J( C )A( £ )r3 (5 »x)dc  • 


(5.16) 
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Here,  J and  A are  the  nucleation  rate  and  the  (effective) 
nozzle  cross-sectional  area,  respectively,  with  both 
as  functions  of  the  position  along  the  nozzle.  The 
quantity  r(£,x)  is  the  final  radius  at  x of  a droplet 
nucleated  at  £ , and  it  is  given  formally  by: 

.x 

* t 

x=e 


r(C,x)  = rU,x)|  + j*  dx  . (5.17) 


The  differential  3r/3x  is  related  to  the  droplet  growth  rate 
by: 

dr/dt  = u(3r/3x)  , (5.18) 

for  steady  flow.  The  first  term  on  the  right  hand  side 
of  Eq.  (5.17)  may  be  identified  as  the  initial  radius 
of  the  droplet  at  its  formation,  i.e.; 

r(C  ,X> I = r (C)  . (5.19) 

'x=c 


The  integrations  in  Eqs . (5.16)  and  (5.17)  are  to  be 
performed  along  the  particle  path  which  in  our  case  is 
equivalent  to  the  nozzle  axis.  The  physical  foundation 
of  Eqs.  (5.16)  and  (5.17)  have  been  explained  before. 


Nucleation  Rate  Equations : 


Despite  the  current  controversy  over  the  physical 


model  of  homogeneous  nucleation  theory,  the  so-called 
"classical  theory" , developed  by  Becker  and  During 
(1935),  Volmer  (1939),  Frenkel  (1946),  et  al . , is 
found  to  give  remarkable  qualitative  agreement  with  many 
experiments  as  will  be  seen  in  Section  6.  For  this  reason  we 
chose  the  classical  theory  as  the  basis  for  incorpora- 
tion in  calculations  performed  with  Eq.  (5.16).  A con- 
venient form  for  this  steady  state  nucleation  rate  is 
given  by: 


■x-JiS-fcK 


n*£n(p  / p ) 


J.,.  (5.20) 


where 


n*  = (4iT/3)p(N./y)(r*)3  , 

C V 


(5.22) 


r*  = (2ooo^Jv)/CRTpcAn(Pv/P»):i 


(5.23) 


In  these  equations,  the  quantities  with  an  asterisk  per- 
tain to  clusters  of  the  critical  size  given  by  Eq.  (5.23), 

■J* 

the  well-known  Gibbs-Thomson  equation.  It  can  be  shown  that 
clusters  of  this  critical  size  are  in  unstable  equilibrium 
with  the  vapor  and  this  size  must  be  achieved 
to  ensure  the  probability  of  further  unlimited  growth  of 
a cluster.  In  cgs  units  Eq.  (5.20)  reduces  to 
T ! 

Confusion  with  respect  to  the  asterisk  assigned  to  states 
at  the  nozzle  throat,  etc.  must  be  avoided. 
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J , = 5 . 4x10 1 9 
cl 


(?) 


(a~Mv) 


1 / 2 


exp 


( n*fn(p  v/Poo)  j 

{'  2 / 


(5.24) 


The  liquid  drop  model  of  the  classical  homogeneous 
nucleation  theory,  i.e.  the  so-called  capillarity 
assumption,  has  been  under  question  for  some  time. 

The  work  of  Kuhrt  (1952a,  1952b),  Lothe  and  Pound  (1962), 

Dunning  (1965),  and  Reiss  (1970),  represents  various 
models  to  incorporate  the  contributions  from  statistical  mech- 
anics considerations  to  the  energy  of  formation  of  the  cluster  and, 
hence,  to  the  rate  of  homogeneous  nucleation.  The  dif- 
ferent ideas  have  resulted  in  rates  which  may  vary  by 
large  factors.  We  shall  not  discuss  the  relative  merits 
of  the  competing  theories.  Rather,  interested  readers  are  re- 
ferred to  the  various  viewpoints  in  the  readily  available  review 
articles  of  Dunning  (1969)  and  Lothe  and  Pound  (1969)  in  the  first 
volume  of  this  series  of  monographs,  and  those  by  Wegener 
and  Parlange  (1970),  and  Reiss  (1970). 

Here,  we  shall  simply  write  down  the  final 
results  of  the  formalism  in  a convenient  general  form 
which  may  be  used  to  compute  the  nucleation  rate  according 
to  any  of  the  models  presented.  Assuming  that  partition 
functions  for  tranlation  and  rotation  of  the  clusters  can 
be  computed  in  the  standard  manner,  we  write: 

cl 


J = r J 


(5.25) 
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where 


r = 


f / 2Trn*m  kT 


(5.26) 


and 


(5.27) 


is  the  moment  of  inertia  of  a spherical  critical  cluster. 

The  "replacement  factor",  , is  introduced  in  the  above 

equation  (Lothe  and  Pound  1966)  to  account  properly  for  the 
changes  in 

energy  in'  the  internal  degrees  of  freedom  of  the  conden- 

* + 

sate  when  a cluster  is  formed  from  the  bulk.  Primarily  solu 

tions  for  the  models  for  q*  are  the  subject  of  the  cur- 

nrep  J 

rent  debate  and  this  factor  is  thus  left  open  in  our  method 
of  calculation. 

On  the  other  hand,  the  validity  of  the  concept  of  a 
macroscopic  surface  tension  for  the  small  clusters  has 
been  questioned.  Since  the  cube  of  the  surface  tension 
of  the  cluster  enters  the  exponent  of  Eq.  (5.20)  via  Eqs . 
(5.22)  and  (5.23),  its  importance  is  self-evident.  Theo- 
retical treatments  of  surface  tension  of  small  clusters 
(Tolman  1949,  Kirkwood  and  Buff  1949,  and  Oriani  and 
Sundquist  1963,  Benson  and  Shuttleworth  1951)  lead  to 

^ The  coefficient  T in  Eq.  (5.25)  will  later  appear  again 
in  Eq.  (6.1)  as  an  empirical  factor  to  match  experiment 
and  theory. 
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sharply  contrasting  results.  However,  we  consider  it 

to  be  consistent  with  the  definition  of  q to  include 

nrep 

in  q the  possible  effects  of  variation  of  surface  tension 
^rep 

with  size.  Thus  we  set  a = am,  and  apply  the  bulk  value  of  the 
surface  tension  in  the  rate  equations  relegating  all  uncertain- 
ties to  the  coefficient  T. 

Thus  we  consider  T as  the  theoretical 
correction  factor  for  the  classical  homogeneous  nucle- 

ation  rate  to  account  for  the  contribution  to  the  energy 
of  formation  of  a cluster  from  statistical  mechanical  as 
well  as  thermodynamic  sources.  Volmer's  work,  etc., 
implies,  of  course,  T = 1. 

In  addition  to  the  steady  state  assumption  Eqs. 

(5.20)  and  (5.25)  are  both  based  on  the 

assumption  that  the  sub-critical  clusters  are  in  thermal 
equilibrium  with  the  ambient  gas  during  their  growth 
towards  the  critical  size.  We  recall  that  this  assump- 
tion is  valid  only  in  the  case  when  a large  quantity  of 
inert  carrier  gas  is  present  to  facilitate  the  dissi- 
pation of  latent  heat  of  condensation  from  the  cluster. 

In  reality,  however,  the  finite  heat  transfer  rate  ought 
to  be  taken  into  consideration  in  nucleation  kinetics 
calculations  to  obtain  a nonisothermal  nucleation  rate,  J 
This  effect  was  first  noted  by  Kantrowitz  (1951),  and  a 
more  detailed  treatment  was  given  by  Feder  et  al . (1966) 
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who  modified  the  classical  kinetic  treatment  of  homo- 
geneous nucleation  to  account  for  the  temperature  dif- 
ferences between  the  cluster  and  ambient  gas  due  to 
finite  heat  transfer  rates.  Their  result  may  be  sum- 


marized as: 


where 


J . = [b2 / (b2  + d2)]  J.  , 

noniso  iso  ’ 


y +1  (1-y  ) m y . + 1 
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(5.28) 


(5.29) 


(5.30) 


and  J.  is  the  '’isothermal"  nucleation  rate.  For 
iso 

J.  we  may  choose  the  classical  rate,  Eq.  (5.20)  or  the 
iso  J 5 ^ 

rate  derived  from  the  statistical  mechanics  treatment, 
Eq.  (5.25). 


Droplet  Growth  Equations : 


After  a cluster  of  the  critical  size  given  by  Eq.  (5.23) 
has  been  formed  by  nucleation,  its  further  growth  by  vapor 
condensation  on  its  surface  is  possible  until  thermo- 
dynamic equilibrium  is  attained.  We  have  mentioned  that 
the  early  stage  of  droplet  growth,  which  is  crucial  in  the 
determination  of  the  onset  of  condensation,  takes  place 


under  free  molecule  flow  conditions.  Resistance 

to  droplet  growth  in  FMF  arises  solely  at  the  gas- 
droplet  interface  and  the  transfer  processes  across 
this  interface  are  characterized  by  the  mass , momentum 
and  thermal  accommodation  coefficients.  This  situation  con- 
trasts with  droplet  growth  under  continuum  flow  con- 
ditions (Kn<<l),  where  the  primary  barrier  to  growth  is 
given  by  vapor  diffusion  and  heat  conduction.* 

The  saturation  vapor  pressure  p of  a droplet  of 

r | 
a pure  liquid  with  radius  r is  given  by  solving  the 

Gibbs -Thomson  equation  Eq.  (5.23)  for  the  vapor  pressure: 

Pr(T)  = p^CTjexp  | 2oyv/ (RTpcr) | . (5.31) 

This  and  the  other  characteristic  properties  of 

droplet  growth  are  indicated  schematically  in  Fig.  5.1. 

Thus,  a droplet  of  size  r is  in  equilibrium  with  the 

surrounding  vapor  if  the  temperature  of  the  droplet  and 

vapor  are  equal  and  if  the  partial  pressure  of  the  vapor 

pv  is  equal  to  p^.  However,  this  equilibrium  is  unstable. 

A larger  droplet  will  tend  to  grow  whereas  a smaller  one 

* In  accordance  with  the  assumption  of  no  slip  between 
the  droplets  and  the  gas  stream,  we  need  to  consider  only 
droplet  growth  in  a quiescent  medium,  i.e.,  bulk  convective 
motion  of  the  surrounding  gas  is  absent. 
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will  evaporate  in  the  same  surroundings. 

The  condensation  of  vapor  molecules  on  the  surface 
of  a droplet  results  in  the  release  of  latent  heat  of 
condensation  which  is  then  dissipated  from  the  droplet.* 
Therefore,  a droplet  growing  in  a supersaturated  environ- 
ment will  generally  have  a surface  temperature  higher 
than  that  of  the  surrounding  medium,  but  lower  than  the 
saturation  temperature  of  the  droplet  corresponding  to 
the  ambient  vapor  pressure.  The  vapor  pres- 
sure exerted  by  the  droplet  will  then  lie  between  the 
free  stream  vapor  pressure  p , and  the  saturation  vapor 
pressure,  p^CT),  at  temperature  T over  a spherical  sur- 
face of  radius  r,  as  given  by  Eq.  (5.31).  These  limits 
are  also  depicted  in  Fig.  5.1.  The  overall  mass  exchange 
rate  on  a droplet  arises  from  the  difference  between  the 
rates  at  which  vapor  molecules  arrive  (condensation)  and 
leave  (evaporation)  the  surface. 

Droplet  growth  behavior  under  continuum  flow  con- 
ditions has  been  studied  by  many  investigators , and  it 
has  been  reviewed  by  Hidy  and  Brock  (1970).  Of 

these  results,  a convenient  formula  for  droplet  growth 

* We  recall  that,  for  many  liquids  at  their  normal  boil- 
ing points,  the  latent  heat  of  condensation  per  molecule 
is  about  12  kT  (Trouton's  rule). 
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rate  in  a quiescent,  continuous  medium  obtained  by 
Gyarmathy  (1963)  is  attractive  in  our  situation.  He 
solved  simultaneously  the  equations  of  vapor  diffusion 
and  heat  conduction  between  the  droplet  and  the  sur- 
rounding medium  and  expressed  the  droplet  growth  rate 
as  a function  of  the  droplet  radius , the  temperature , 
pressure  and  the  composition  of  the  surrounding  gas. 

The  fact  that  the  surface  temperature  of  the  droplet 
does  not  appear  explicitly  in  Gyarmathy' s formula  for 
droplet  growth  in  a continuum  makes  his  results  more  use- 
ful than  other  results,  since  the  surface  tempera- 
ture is  unknown  a priori. 

The  droplet  growth  rate  in  FMF  may  be  expressed 
as 

^ £D(pv,T)-DeCTd,r)J  , (5.32) 

where 

D(p  ,T)  = p ( 2irm  kT)"*5  , (5.33) 

is  the  impingement  rate  of  vapor  molecules  from  kinetic 
theory,  De(Td,r)  is  the  evaporation  rate  of  a droplet  at 
temperature  T^  and  radius  r,and  a is  the  mass  accommo- 
dation coefficient.  Under  the  assumptions  underlying 
the  principle  of  detailed  balancing,  this  evaporation 


rate  may  be  taken  (Knudsen  1915,  Paul  1962)  to  be 


equal  to  the  impingement  rate  on  the  droplet  by  its 


saturated  vapor  at  the  same  temperature. 


Thus  : 


De(Td,r)  = D [pr(Td),  Td  ] = 


*r'V 


LI  . 

c d 


, (5.34) 


where  pr(Td)  is  again  given  by  Eq.  (5.31).  Alterna- 
tively, we  may  define  an  effective  condensation  coef- 
ficient that  formally  accounts  for  evaporation  by 


a'  = a [l-De(Td,r)/D(pv,T)J  , 


(5.35) 


and  rewrite  Eq . (5.32)  as: 


dr/dt  = a'(m  /p  )D(p  ,T) 
c c v 


(5.36) 


Since  the  temperature  of  the  drop*  is  unknown,  it 
must  be  determined  by  energy  considerations.  Hill  (1966) 
derived  the  following  equation  for  the  energy  balance  of 
this  droplet; 


ctK  k(DT-De  T,)/rn  + (l-ct)K  kD£  (T-T,)/m  + D.  K.  . (T-T  . )/m. 

v dc  vv  dc  il  l d a 


«h  (T , )aD(  1-De/D)  , 

c d ’ 


(5.37) 


* We  assume  the  temperature  inside  these  droplets  to  be 
uniform.  It  can  be  shown  that  the  temperature  equali- 
zation time  inside  the  small  droplets  under  consideration 
(r<1000A)  is  extremely  short,  of  the  order  of  0.1  usee 
or  shorter. 
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where  K = ( y+1 ) / [ 2 ( y-1 ) ] so  that  the  product  KkT  denotes 
the  average  energy  of  impinging  molecules . The  term 
hc(Tc})  = Eyv/(yv-1)  ]kTd/mc-L(Td)  is  the  enthalpy  of  the 
condensate  at  temperature  Td . The  first  term  on  the 
left  hand  side  accounts  for  the  difference  in  energy 
carried  by  the  condensing  and  evaporating  molecules. 

By  definition  of  the  thermal  accommodation  coefficient 
£,  the  second  and  third  terms  denote  the  energy  trans- 
ferred by  the  reflection  of  vapor  and  inert  carrier  gas 
molecules,  respectively,  at  the  droplet  surface.  The 

I 

term  on  the  right  hand  side  is  the  rate  of  energy 
accumulation  of  the  drop  neglecting  the 
relatively  small  contributions  from  temperature  change 
and  surface  energy  of  the  drop  (Hill  1966). 

Together  with  Eq.  (5.31)  for  the  equilibrium  vapor 
pressure  over  a curved  surface,  Eqs.  (5.35),  (5.36)  and 
(5.37)  may  be  solved  simultaneously  for  the  three  unknowns 
dr/dT,  Td,and  p^(Td).  Owing  to  the  complexity  of  these 
equations,  an  iteration  procedure  is  necessary  to  obtain 
the  growth  rate  dr/dt,  (e.g.,  Hill  1966,  Jaeger  et  al. 

1969,  Dawson  et  al . 1969). 

Fortunately,  for  most  cases  of  practical  interest, 
such  an  exact  solution  is  not  necessary  and  useful  approxi- 
mate solutions  have  been  obtained  by  Wu  (1972).  For 
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example,  when  a large  quantity  of  inert  carrier  gas 
is  present,  the  droplet  remains  practically  in  thermal 
equilibrium.  Typically  in  moist  air,  a water  molecule 
undergoes  roughly  a thousand  collisions  with  inert  gas 
molecules  for  each  encounter  with  another  water  molecule  which  may 
lead  to  condensation.  Hence,  = T,  the  condition  of 
the  droplet  is  represented  by  point  2 in  Fig.  5.1,  and 
we  have : 

/ 

(5.38) 


a'  = a[l-pr(Ty'pvj 


Next,  when  the  droplet  and  the  ambient  gas  are  nearly 
in  thermal  equilibrium  i.e.  its  state  is  near  point  2 
in  Fig.  5.1,  an  approximate  expression  for  a'  has  been 
obtained : 


a Cl  - Pr<T)/pv] 

a'  = p (T)  ~Lm  Lm  /kT  -l/ 2 

, j a *r  c c 

1 K p kT 
v *v 


(5.39) 


ci-d-a)  C1-5V)  J + Ci  Y£i7M3/2 

Kv  Pv\mi/ 


Equation  (5.39)  is  valid  only  when  | T^-T | /T<<1  , however, 
under  practical  conditions  this  requirement  is  satisfied  when 


kTp 


m L 
c 


p m L 
*v  c 


kT 


>>1 


(5.40) 


which  is  normally  the  case, 
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On  the  other  hand,  for  a high  latent  heat  of  con- 
densation of  the  vapor,  or  a low  concentration  of  inert 
carrier  gas , the  droplet  temperature  approaches  the 
saturation  temperature  for  the  droplet,  indicated  by 
point  3 in  Fig.  5.1.  In  the  limit  when 

= Tg  (r,pv),  the  saturation  temperature  at  pv  over 
a droplet  of  radius  r,  we  have,: 

a'  = o(l-/T7TT)  • (5.41) 

s 

For  droplet  temperatures  close  to  the  saturation  tem- 
perature , | Td-Tg  | <<Ts,  we  have: 


a 


a'  = 


L-Vi-  r i- 

s L 
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(5.42) 


[l-(l-a)  (1-CV) ] + 5 
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Equations  (5.38),  (5.39),  (5.41)  and  (5.42)  cover  a sub- 
stantial part  of  the  entire  temperature  range  of  a 
growing  droplet.  In  addition,  the  "near  equilibrium" 
and  "near  saturation"  approximations,  Eqs . (5.39)  and 
(5.42)  may  also  be  used  to  compute  the  evaporation  rate 
of  a droplet  in  an  unsaturated  atmosphere,  resulting  in  a 
negative  a'.  These  equations  are  explicit  ex- 

pressions for  a*  in  terms  of  measurable  ambient 


■ 
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conditions  only,  therefore  they  offer  an  attractive 
alternative  to  solving  the  mass  and  heat  transfer 
equations  simultaneously. 

Notice  that  the  terms  entering 
Eq.  (5.16)  and  the  nucleation  and  droplet  growth  rate 
equations  are  generally  functions  of  the  thermodynamic 
state  of  the  two-phase  system  under  consideration.  The 
state  of  the  system  is  determined  by  the  gasdynamics 
of  the  flow,  which,  in  turn,  is  influenced  by  the  con- 
densation process.  It  is  clear  that  the 
kinetics  of  phase  change  is  coupled  with  the  gasdynamics 
of  the  flow  process.  Equations  (5.16)  and  (5.17)  must 
therefore  be  integrated  simultaneously  with  the  proper 
gasdynamic  equations  of  motion  and  the  mass  fraction  of 
the  condensate  will  then  emerge  as  part  of  the  complete 
solution  . * 

6 EXPERIMENT  AND  THEORY 

We  can  now  discuss  typical  results  on  homogeneous 
nucleation  as  obtained  by  gasdynamic  methods.  At 
first,  general  examples  will 

* Computer  programs  have  been  written  at  our  laboratory 
(Wu  1974)  for  this  problem  and  several  variants.  They 
are  available  upon  request. 


be  given  on  the  functional  relationships  between 
experiment  and  theory.  Next,  a data  collection  will 
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be  provided  with  many  results  newly  evaluated  by  iden- 
tical methods  of  data  reduction.  While  we  rely  largely 
on  information  from  our  laboratory  in  the  first  part, 
the  tabular  results  on  several  substances  which  have 
not  been  published  previously  will  include  experiments 
from  other  sources  provided  that  the  "raw  data"  could 
be  retrieved  from  the  publications.  It  is  believed  to 
be  imperative  to  treat  experimental  findings  uniformly 
because  of  the  great  sensitivity  of  various  parameters 
in  nucleation  theory  to  the  choice  of  various  properties. 
Finally,  all  comparisons  have  been  screened  to  comply 
with  the  assumptions  implicit  in  the  previous  sections 
including  the  following  criteria. 

One-dimensional  steady  or  unsteady  flow  prevails. 
Viscous  effects  and  turbulence  in  the  free  stream 
are  unimportant . 

The  flow  parameters  are  known  throughout  the  conden- 
sation region  (or  at  least  at  the  onset  of  conden- 
sation) . 

The  heat  addition  by  condensation  does  not  lead  to 
thermal  choking,  i.e.  no  shock  waves  appear. 


& 
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Continuum  flow  prevails.* 

Cooling  rate  and  time  scales  are  sufficiently  mod- 
erate for  the  steady-state  assumption  of  the  classi- 
cal nucleation  theory  to  apply.* 

The  presence  of  aerosols  is  unimportant. 

Vapor  impurities  are  of  sufficiently  low  level 
that  binary  nucleation  is  unimportant, 
a)  General  Results. 

Properties  at  the  onset  of  condensation.  Flow  experi- 
ments provide  raw  data  such  as  those  of  Fig.  4.1  directly 
yielding  the  onset  conditions  as  shown  in  Fig.  4.4.  An  under- 
standing of  the  state  parameters  through  the  condensation 
zone  (Fig.  4.2)  can  moreover  be  derived  by  solving  the  diabatic 
flow  equations  with  heat  addition  as  described  in  Section  3b.** 
Since  the  heat  addition  is  due  to  condensation,  Eq.  (3.19) 
can  be  applied  to  the  data  of  Fig.  4.2.  Thus,  the  conden- 
sate mass  fraction  as  a function  of  distance,  (or  time, 

(see  Eq.  (2.32)),  can  be  determined  directly  as  shown  in  Fig. 

6.1.  The  shaded  area  describes  the  uncertainty  of  g(x) 
resulting  from  experimental  uncertainties.  It  is  this 
remarkable  resolution  of  the  rate  of  condensation  beyond 
the  critical  supersaturation  that:  sets  the  nozzle  method 
apart  from  the  cloud  chamber.  However,  g(x)  cannot  yet  be 


* This  assumption  again  rules  out  the  inclusion  of  cluster- 
ing in  high  intensity  molecular  nozzle  beams  e.g.  Hagena 
(1974). 

**  The  variations  in  molecular  weight,  etc.  of  the  mixture 
as  a result  of  change  in  mixture  composition  may  be  incor- 
porated in  the  equation  without  difficulty. 
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computed  independently  based  on  first  principles  by  the 
inclusion  of  nucleation  and  growth  equations . To  apply  the 
equations  developed  in  Section  5 for  comparison  it  is 
therefore  useful  to  introduce  an  empirical  adjustment 


factor,  T,  such  that 


r = J /J.v  , 
exp  th  * 


(6.1) 


where,  symbolizes  a theoretical  expression  for  the 

nucleation  rate  such  as  Volmer's  Eq.  (5.24)  or  some  other 
expression  for  the  nucleation  rate  such  as  Eq.(5.28).  For  a 
given  experimental  environment  and  nozzle 
geometry,  a complete  solution  of  the  equations  can  be 
repeated  for  an  estimated  range  of  values  of  r in  Eq. 

(6.1).  That  value  of  T can  then  be  chosen  for  which 
the  strong  increase  in  computed  nucleation  rate  coin- 
cides with  the  observed  onset  of  condensation.  This 
location  in  space  (nozzle)  or  time  (shock  tube),  can 
be  more  clearly  defined  moreover  as  a match  of  the 
observed  and  the  computed  condensate  mass  fraction  e.g. 
at  the  location  of  the  departure  of  the  pressure  from 
the  isentrope  in  Fig.  4.1. 

We  recall  that  the  symbol,  T,  was  previously 
introduced  in  Eqs.  (5.25)  and  (5.26)  as  a ratio  of  parti- 
tion functions.  The  empirically  determined 

factor,  T,  of  Eq.  (6.1)  as  obtained  from  a complete 
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solution  of  nucleation  and  growth  equations  in  compari- 
son with  experiment  can  now  be  interpreted  as  an  experi- 
mental solution  for  Eq.  (5.26).  This  proposal  will  be 
elucidated  in  the  following  Section  6b.  At  the  onset 
of  condensation  here  defined  as  the  location  in  Fig.  6.1  where 
g = 0.7x10  3ft,  we  note  F = 10 5 for  a = in  Eqs . (5.20) 
and  (5.25).  For  ethanol,  the  temperature  at  the  onset 
of  condensation  lies  above  the  triple  point  permitting 
the  use  of  tabulated  surface  tension  values  for  the  bulk 
liquid.  In  this  fashion  different  versions  of  the  theory 
of  nucleation  mav  readily  be  compared  with  experiment. 

The  well-known  "go"  versus  "no  go"  feature  of  nucle- 
ation is  again  demonstrated  in  Fig.  6.1  since  the  nucle- 
ation rate  plotted  quickly  attains  a maximum  of  about 
J~1019  cm” 3 sec-1  at  the  measured  onset  of  condensation. 

About  40  ysec  before  the  vapor  sample  reaches  this 
location,  the  rate  of  nucleation  was  as  low  as  J = 1 per 
earth  volume  per  age  of  the  earth!  While  the  cloud  chamber 

experiment  may  at  best  establish  the  critical  supersatu- 
ration of  the  onset  in  an  undisturbed  small  central  volume 

* This  relatively  higher  value  of  the  condensate  mass  frac- 
tion than  those  used  later  for  water  vapor  is  required 
since  the  latent  heat  of  ethanol  is  lower  and  consequently 
more  condensate  is  needed  to  measurably  disturb  the  isen- 
tropic  flow  by  heat  addition.  In  fact, if  a diagnostic 
method  such  as  light  scattering  permits  the  location  of  the 
onset  of  condensation  experiments  for  vanishing  heat  addi- 
tion can  be  performed  permitting  the  application  of  the 
isentropic  properties . 


in  the  chamber,  the  "sensitive  volume"  (Kassner  and 
Schmitt  1966),  the  gasdynamic  evaluation  permits  a 
detailed  resolution  of  this  rapid  process. 

Properties  of  the  condensation  zone.  Once  the 
measured  onset  conditions  have  been  matched  by  theory , 
a further  complete  calculation  of  nucleation  and  growth 
yielding  the  accretion  of  the  condensate  mass  fraction 
beyond  the  onset  location  can  be  computed  from  data  such 
as  those  of  Fig.  4.2.  These  results  for  g(x)  agree  well 
with  the  measurement  as  shown  in  Fig.  6.1.  The  calcu- 
lation of  droplet  growth  with  a = 0.7  in  Eqs.  (5.3“6)  and 
(5.39)gives  excellent  results.*  In  this  instance  nucle- 
ation and  growth  follow  each  other,  each  taking  about 
20  ysec.  We  note  the  fast  build-up  of  the  number  of  parti- 
cles and  the  rapid  drop  of  the  nucleation  rate  resulting 
from  the  decline  in  saturation  ratio  in  the  condensation  zone 
caused  by  vapor  depletion.  The  number  density  of  the  parti- 
cles again  drops  slightly  due  to  the  ensuing  expansion 
while  coagulation  can  be  shown  to  be  negligible  in  the 
condensation  zone  (Wegener  and  Stein  1969).  Finally,  at 
the  end  of  the  condensation  zone,  thermodynamic  equilibrium 


* The  choice  of  a = 0.7  is  reasonable  in  the  light  of 
recent  work  on  this  coefficient  (Mills  and  Seban  1967, 
Lednovich  1974). 


is  again  established  when  nearly  all  the  ethanol  has  conden- 
sed, i.e.,  in  mg.  In  fact  the  partial  pressure  of 
ethanol  vapor  required  to  maintain  Pv/pro  = 1 is  negligible 
at  these  low  temperatures. 

With  the  same  methods, we  can  now  interpret  light 
scattering  results  in  greater  depth  than  shown  in  Fig. 

4.3.  Computed  distribution  functions  for  cluster 
populations  emerge  as  by-products  of  these  calculations 
as  plotted  in  Fig.  6.2.  Most  of  the  droplets  are  roughly 
of  the  same  size  reflecting  the  fact  that  most  of  them 
originated  on  nuclei  of  critical  size  formed  at  the 
location  of  the  peak  value  of  the  nucleation  rate 
slightly  beyond  the  onset  location.  Similar  results 
are  available  for  H2O  in  air  (Wegener  and  Stein  1969), 
and  pure  steam  (Barschdorff  et  al.  1974).  In  all  situations 
the  computed  and  measured  drop  sizes  agree  well. 

Cooling  rate.  At  the  high  cooling  rates  in  rapid 
expansions  the  question  arises  if  the  steady-state  nucle- 
ation rate  can  be  established  in  the  short  time  available. 
These  induction  times  have  been  variously  estimated  as 
reviewed  e.g.  by  Courtney  (1962),  and  Andres  and  Boudart 
(1963  ).  Ordinarily  , such  time  scales  are  estimated  to  be 
of  the  order  of  a few  microseconds,  however,  Courtney 
(1962)  suggests  larger  values  at  the  low  temperatures  in 
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question.  At  any  rate  with  gasdynamic  techniques,  again 
with  the  striking  exception  of  the  high  intensity  molecu- 
lar beams , the  time  of  persistence  of  the  supersaturated 
state  and  the  time  elapsing  in  the  actual  condensation 
zone,  e.g.  Fig.  6.1,  appear  to  be  significantly  above 
the  induction  time.*  We  can  demonstrate  this  fact  with- 
out resort  to  basic  theory  by  calculations  on  the  effect 
of  cooling  rate  for  condensation  in  two  nozzles  of  dif- 
ferent cooling  rates.  Figure  4 . la  shows  an  experiment  in 
a nozzle  with  a cooling  rate  computed  from  Eq.  (2.34)  of 
dT/dt  = - 0.7  C/ysec.  With  an  empirical  value  of  log1Q  T 

= 4.7  and  a = 0.6  in  Eqs . (5.25).  (5.36),  and  (5.39) 
respectively,  together  with  the  equations  of  motion  with 
heat  addition  we  obtain  the  calculated  results  for  pressure, 
also  shown  in  Fig.  4.1a.  This  calculation  again  gives 
excellent  agreement  between  experiment  and  theory.  Next 
we  apply  the  same  values  of  T and  a without  any  further 
adjustment  to  a second  complete  calculation  for  the  con- 
ditions of  the  experiment  shown  in  Fig.  4.1b.  Here  we 

deal  with  a nozzle  with  a cooling  rate  of  dT/dt  = - 2.8  C/ysec, 

* A DamkBhler  parameter  giving  the  ratio  of  a character- 
istic flow  time  to  the  induction  time  or  the  time  elapsed 
in  the  condensation  zone  offers  itself.  Such  a discussion 
is  of  particular  interest  in  the  nozzle  flow  treated  as  a 
nonequilibrium  flow  with  a single  relaxing  mode  without 
regard  to  the  physical  mechanism  of  condensation  (Wegener 


1975)  . 
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and  again  excellent  agreement  with  the  experimental 
result  is  achieved.  Although  the  cooling  rates  differ  only 
by  a factor  of  four  between  the  two  nozzles,  it  appears  that 
the  theory  yields  satisfactory  results  in  this  range 
of  expansion  rates.  If  the  steady-state  assumption 
were  not  valid  in  these  rapid  expansions , such  agree- 
ment would  be  quite  unlikely. 

Approach  to  the  critical  point . Nucleation  theory 
implies  (Volmer  1939)  that  an  initial  step  towards  the 
build-up  of  critical  nuclei  involves  the  formation  of 
a dimer  by 

A+A+MtA„+M,  (6.2) 

z 

where  M represents  a third  body.  The  start  of  nucle- 
ation is  therefore  tied  to  tertiary  collisions.  It  is 
interesting  to  note  that  in  fact  in  the  onset  region  of 
typical  nozzle  experiments  such  as  shown  in  Figs.  4.1 
or  6.1,  triple  collisions  are  about  one  order  of  magni- 
tude more  frequent  then  the  nucleation  rate  in  this  range. 
With  the  collision  frequency  depending  strongly  on  pres- 
sure, it  is  expected  that  the  supersaturation  decreases 
at  increased  pressure,  in  particular  if  the  critical 
point  is  approached.  Moreover,  the  surface  tension 
decreases  to  vanish  at  the  critical  temperature.  This 
feature  of  the  theory  is  indeed  borne  out  by  nozzle  and 


shock  tube  experiments  with  steam  as  seen  in  Fig.  6.3. 
The  results  of  Gyarmathy  (1973)  at  temperatures 
approaching  the  critical  point  are  particularly  remark- 
able since  they  were  obtained  in  situ  in  a steam  tur- 
bine operating  in  a modern  power  plant  (Karlsruhe, 
Germany)  with  pressures  at  the  condensation  onset  of 
up  to  75  atmospheres.  Additional  evaluation  of  these 
high  pressure  data  is  difficult  since  the  ideal  gas  law 
assumptions  are  inapplicable  for  steam  once  the  pres- 
sures in  the  flow  exceed  about  one  atmosphere. 

Carrier  gas  effects. 

The  isothermal  assumption  implicit  in  the  deriva- 
tion of  Eq.  (5.20)  is  physically  well  realized  by  mixing 
the  vapor  with  an  excess  of  a non-condensing  carrier  gas 
as  found  e.g.,  in  the  moist  air  in  a cloud  chamber.  In 
recent  shock  tube  experiments  on  H2O  condensation  mixed 
in  He,  A,  and  air,  shown  in  Fig.  6.4,  Barschdorff  (1975) 
verified  the  absence  of  carrier  gas  effects.  The  slight 


differences 
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noted  can  be  rationalized  in  view  of  the  different 
values  of  the  ratio  of  the  specific  heats,  y,  of  the 
gas-vapor  mixtures  leading  to  slightly  different  cool- 
ing rates  in  the  centered  expansion  in  the  shock  tube 
driver  chamber  (see  Section  2c).  Recalling  our  dis- 
cussion on  droplet  growth  we  note  that  the  number  of 
collisions  available  exceeds  the  number  required  to 
insure  thermal  equalization  of  cluster  and  surrounding 
gas  mixture  demonstrating  that  the  properties  of  the 
carrier  become  unimportant. 

Non-isothermal  condensation.  The  simplifying  condi- 
tion of  equal  temperature  of  drops  and  vapor  ceases  to 
be  valid  in  pure  vapor  condensation.  Here, the  appli- 
cation of  the  classical,  isothermal  theory,  e.g.  Eq. 
(5.20)  with  a = and  a = 1 to  a complete  solution  of 
the  equations  for  typical  steam  nozzle  flow  in  nozzles 
yields  r “ 10-5,  i.e.  the  experimental  rate  of  conden- 
sation appears  to  be  appreciably  slower  then  that  pre- 
dicted by  Volmer's  equation.  Experiments  of  many 
steam  investigators  compared  with  theory  as  listed  by 
Barschdorff  et  al.  (1972)  prove  this  point  and  new 
calculations  are  given  in  Table  6.2.  On  an  identical 
basis,  the  comparison  of  theory  and  experiment 


(Arthur  1952)  for  pure  nitrogen  yields  T ~ 10-6  also 
seen  in  Table  6.2.  Oddly  enough  until  recently 
(Barschdorff  et  al . 1972), this  fact  appears  to  have  been 
overlooked  excepting  a general  remark  in  the  work  of 
Deych  et  al.  (1969a)  who  pointed  out  that  the  theory 
yields  condensation  rates  for  stehm  "which  are  too  high 
by  several  orders  of  magnitude"  in  comparison  with  ex- 
periments. Usually  this  discrepancy  between  the  theory 
and  the  experiment  was  masked  by  artificially  "slowing 
down"  the  condensation  rate  in  the  theory  by  apply- 
ing low  values  of  the  condensation  coefficient  to  the  growth 
law  (Hill  1966),  or  resorting  to  the  unjustified  introduc- 
tion of  a "friction  factor"  in  the  momentum  equation  (Campbell 
and  Bakhtar  1973).  We  should  in  fact  expect  values 
0.1<a<l  (Mills  and  Seban  1967)  for  the  mass  accommo- 

dation coefficient  and  boundary  layer  effects  can  be  assumed 
to  be  negligible  or  calibrated  via  the  effective  area  con- 
cept (Section  2).  The  apparent  low  values  of  T must 
therefore  be  accounted  for  otherwise. 

Feder  et  al . (1966)  based  on  prior  work  (Kantrowitz  1951) 
treated  the  non-isothermal  build-up  of  clusters  in  pure 
vapors  by  deriving  a new  preexponential  factor  for  the 
nucleation  rate  as  shown  in  Eq.  (5.28).  This  non- 
isothermal  factor  yields  a lower  nucleation  rate  since 
some  of  the  vapor  molecules  that  collide 
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with  the  droplets  do  not  condense  and  they  now  carry 
off  the  thermal  energy  arising  from  the  latent  heat.* 

An  application  of  this  result  to  steam  condensation  is 
demonstrated  in  Fig.  6.5.  At  large  values  of  the  con- 
densation coefficient,  for  the  growth  rate  of  water 
droplets  beyond  their  critical  size , the  rate  of  nucle- 
ation  is  slowed  with  respect  to  the  isothermal  equation 
by  non-isothermal  effects  by  about  three  orders  of 
magnitude.  A proper  combination  of  the  non-isothermal  theory 
and  a reasonable  value  of  a provide  good  estimates  for 
technological  purposes  of  the  onset  of  condensation  in 
steam  nozzles . 

Finally,  additional  features  of  the  experiments 
with  varying  vapor  mass  fractions  (Barschdorff  1975) 
compared  with  the  non-isothermal,  classical  nucleation 
theory  (with  a = am  and  a = 1 in  the  growth  law)  are 
seen  in  Fig.  6.4.  Values  for  Eq.  (6.1)  with  the  non- 
isothermal  theory  are  plotted  as  a function  of  the  ini- 
tial specific  humidity  in  the  range  of  0.004<ajQ<l,  with 
the  previously  discussed  results  from  different  carrier 
gases.  In  the  shock  tube  experiments,  the  cooling  rate 

* This  effect  is,  of  course,  equivalent  to  the  use  of  a 
low  "sticking  coefficient",  however,  the  physical  process 
is  of  a different  origin,  see  Section  5. 
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is  lower  by  about  one  order  of  magnitude  (10 5 vs.  10 6 C/s) 
than  in  the  nozzle  experiments  of  Fig.  6.1.  The  puzzling 
result  of  Fig.  6.4  is  the  fact  of  a variation  of  106<T<1 
between  the  limits  of  low  vapor  concentration  to  pure 
steam.  In  light  of  our  previous  thoughts,  one  might  expect 
T = constant  (although  possibly  T / 1)  irrespective  of  the 
initial  humidity  since  critical  drop  sizes,  expansion  rates, 
initial  conditions,  etc.  are  roughly  equal.  Although  the 
non-isothermal  theory  of  nucleation  predicts  the  proper 
general  trend  as  seen  in  Fig.  6.5,  the  dominant  exponential 
term  in  the  rate  equation  is,  however,  unaltered  in  this 
work.  Here  the  surface  tension  enters  in  the  third  power 
(Parlange  1970).  The  calculations  in  Fig.  6.4  were  per- 
formed for  a a (T)  function  obtained  from  tables  for  the 

OO 

temperature  of  the  gas  flow  as  obtained  for  the  isentropic 
expansion.  However,  the  droplet  is  warmer  then  its  surround- 
ings affecting  its  surface  free  energy.  Since  the  droplet 
temperatures  have  not  been  calculated,  the  trend  noted  in 

Fig.  6.4  remains  quantitatively  unexplained. 
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Small  cluster  properties.  Finally, some  quali- 
tative hints  can  be  provided  on  the  properties  of  small 
clusters  by  making  use  of  the  phenomenological  accuracy  of 
the  classical  theory  of  nucleation.  By  "small"  we  con- 
sider here  the  typical  critical  H2O  drops  or  clusters 
whose  critical  radii  as  computed  with  Eq.  (5.23)  of 
4<r*<12R  for  typical  onset  data  contain  about  9<n*<250 
molecules  (see  Table  6.1)  as  obtained  from  Eqs.  (5.23) 
and  (5.22).  Trends  of  experimental  data  can  be  deduced 
by  assuming  that  the  pre-exponential  in  Eq.  (5.20)  is 
constant  for  a given  expansion  to  get: 

J = const,  exp  £-  const.  (o^/T)3  £n-2 (p/p^) J . (6.3) 

If  we  further  set  a = const,  and  J = const,  at  the  onset 

00 

of  nucleation , we  expect  the  function: 

log(p /pj  = f (t~2/3)  , (6.4) 

to  be  a straight  line.  In  fact  cloud  chamber  results  as 
reproduced  in  Fig.  6.6  verify  this  trend.  We  find 
moreover  log  (p/p^)  = 0 for  T = «>  as  expected  formally. 

In  the  same  figure,  we  show  experimental  data  on  water 
vapor  condensation  in  nozzles  and  shock  tubes  from  a 
variety  of  sources.  For  simplicity,  the  vapor  pres- 
sure was  taken  to  be  that  over  ice  if  T<273  K in 
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order  to  have  a uniform  basis  for  evaluation  of  the 
critical  supersaturation.  Many  other  results  were 
studied  for  H2O  and  other  vapors  in  the  same  manner 
and  in  general  the  trends  are  similar  to  those  of 
Fig.  6.6.  In  line  with  the  empirical  function  of 
Eq.  (1.1)  and  Figs.  1.3  and  1.4,  lower  humidi- 

ties and/or  higher  cooling  rates  respectively  lead  to 
a larger  critical  supersaturation.  Conversely,  the  two 
effects  may  counteract  each  other.  In  general,  it  is 
found  that  given  sets  of  experiments  lie  on  straight 
lines,  however,  the  lines  do  not  point  to  the  origin 
of  the  coordinate  system,  as  formally  required  by  Eq. 

(6.3).  Similar  results  are  noted  for  steam  as  seen  in 
Fig . 6.7. 

A straight  line  is  noted  even  when  the  vapor  be- 
havior departs  drastically  from  the  perfect  gas  behavior 
as  seen  for  the  high  pressure  data  of  Gyarmathy  et  al. 

(1973).  His  results  point  to  the  critical  point  where  we  expect 
log(p/poo)  = 0 at  ^crit  = 374  C.  It  is  inviting  to 
relate  the  slopes  of  these  lines,  their  intercepts,  etc. 
to  additional  features  of  the  theory.  The  results  of 
Figs.  6.6  and  6.7  show  that  the  "constant"  preceding  the 
exponential  on  the  right  hand  side  of  Eq.  (6.3)  is  itself 
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a function  of  the  temperature  for  the  correlation  given. 
Direct  use  of  Eq.  (6.3)  including  a variable  bulk  sur- 
face tension  does  not,  however,  materially  alter  the 
trend  of  the  curves  shown.  It  therefore  appears  that 
the  surface  tension  of  the  clusters  is  not  an  important 
variable  in  this  range  pointing  to  a structure  whose 
properties  depend  little  on  temperature.* 

Binary  Nucleation. 

Finally,  gasdynamic  methods  have  been  applied  just 
like  the  cloud  chamber  to  a study  of  the  binary  nucle- 
ation in  the  ethanol-water  system.  Although  this  sub- 
ject is  outside  our  present  scope,  we  refer  to  nozzle 
experiments  discussed  by  Clumpner  (1970)  and  Wegener 
et  al.  (1972),  and  recent  shock  tube  work  by  Wu  and 
Belle  (1974). 

b)  Experiment  and  Theory  for  the  Critical  Super- 
saturation . 


This  final  section  will  be  concerned  with  quantitative 
comparisons  of  computed  and  observed  critical  supersatu- 
ration. However,  at  first  a remark  on  the  accuracy  of 
such  comparisons  is  required.  We  recall  that 


* Stein  and  Armstrong  (1973)  show  that  e.g.  H2O  clusters 
in  a molecular  beam  exhibit  a diamond  cubic  structure. 
However,  in  these  experiments  at  extreme  cooling  rates 
(see  Section  4c),  the  internal  degrees  of  freedom  and  the 
rotational  mode  may  well  have  been  "frozen"  at  the  state 
in  the  supply.  Thus  a comparison  with  the  conditions 
discussed  here  is  not  possible. 


depending  on  the  latent  heat  a certain  value  of 
the  condensate  mass  fraction  is  detectable  by  gasdynamic 
methods.  Consequently  values  of  T in  the  theory,  see 
Eq.  (5.26)  or  empirically,  see  Eq.  (6.1)  are  matched  at 
this  point  of  departure  from  isentropic  flow.  In  con- 
duction with  an  error  analysis  of  all  experimental 
uncertainties,  it  is  possible  to  determine  the  resulting 
error  in  the  measured  nucleation  rate.  Such  an  analysis 

s 

was  made  for  ethanol  condensation  in  a steady  nozzle 
flow  (Clumpner  1970).  Here  the  detectable  minimum  mass 

V 

fraction  is  larger  than  that  for  water  owing  to  the 
lower  latent  heat  and  therefore  the  errors  may  be  larger 
leading  to  a conservative  result.  Typically,  a pressure 
deviation  from  the  isentrope  of  about  four  times  the 
experimental  error  leads  to  a detection  of  a mass  frac- 
tion oj*g£l0-3.  The  result  of  the 'analysis  reveals  that 
at  the  onset  location  of  condensation  in  the  nozzle  in 
a typical  ethanol/air 


experiment  (e.g.  Fig.  6.1),  the  nucleation  rate  can  be 
determined  within  about  ±20%.  Somewhat  better  values  are  to  be 
expected  for  water  where  the  higher  latent  heat  permits 
the  detection  of  smaller  amounts  of  condensate.  However, 


roughly 
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corresponds  to  the  ideas  initially  advanced  by  Kuhrt 
(1952),  and  by  Dunning  (1965)  leading  to  "low"  values  of  r 
in  contrast  to  the  work  of  Lothe  and  Pound  (1966). 

Plotting  the  main  terms  in  the 

nucleation  rate  equation  during  the  course  of  the  entire 
condensation  zone  of  Fig.  6.1  we  find  Fig.  6.9.  Only 
the  classical  exponential  term, here  determined  for 
a = a , varies  dramatically  during  the  expansion.  The 
factor  T of  Eq.  (5.26)  computed  with  the  Dunning  model 
for  q*rep  agrees  well  with  the  experiment.  Such 
graphs  permit  detailed  comparisons  with  the  theory  not 
offered  by  other  experimental  techniques.  In  particular, 
work  of  this  kind  ought  to  be  extended  to  other  substances 
with  other  properties . 

Tables  6.2  and  6.3  finally  give  results  from  theory 
and  experiment  matched  at  the  critical  supersaturation. 

The  non-isothermal  theory  and  the  isothermal  expression 
are  compared  for  pure  vapor  expansions  in  Table  6.2.  In 
Table  6. 3, only  the  isothermal  expressions  are  used  since 
the  condensing  vapors  were  expanded  in  an  excess  of  car- 


rier gas . 


However,  once  condensation  was  observed  at 


* We  recall  the  previous  observations  on  the  small  effect 
of  the  nature  of  the  carrier  gas  on  the  onset  of  condensa- 
tion. 
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temperatures  below  the  triple  point,  vapor  pressures 
with  respect  to  the  liquid  (estimated  for  supercooled 
states),  or  solid  phase  respectively  were  computed 
leading  to  distinctly  different  answers.  Table  6.1  aids 
in  relating  n*  to  the  values  of  r*  of  Tables  6.2  and 
6.3.  The  tables  are  self-explanatory  and  they  permit 
varying  interpretations.  However,  the  results  shown 
are  believed  to  be  the  first  comparative  data  collec- 
tion uniformly  based  on  identical  criteria  as  discussed 
in  Section  5. 

A few  observations  on  the  results  may  be  in  order. 

For  one,  the  critical  radii  of  the  clusters  in  the 
range  of  4<r*<llX  are  appreciably  smaller  then  those 
observed  in  cloud  chambers.  In  the  latter  experiments 
lower  supersaturation  is  experienced  due  to  the  lower 
cooling  rate  in  line  with  the  empirical  relation  Eq. 

(1.1).  The  small  size  of  the  clusters  (see  Table  6.1) 
may  invalidate  all  comparisons  with  the  classical 
theory  and  its  variants.  Moreover, for  these  small 
clusters  the  distinctions  between  the  liquid  and  solid 
phase  and  corresponding  property  selections  become 
blurred  as  noted  before.  In  fact  it  is  impossible  at  present 

to  distinguish  a priori  if  liquid  or  solid  properties  apply. 

A case  in  point  is  benzene  as  noted  in  lines  1 and  2 of 
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A 
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Table  6.3.  Applying  solid  properties  at  60  to  80  C 
below  the  triple  point  leads  to  a supersaturation 
similar  to  that  of  water.  In  sharp  contrast,  the 
assumption  that  the  condensate  is  a supercooled  liquid 
yields  a nucleation  rate  of  about  10*  times  that  of 
the  classical  rate,  i.e.  r-»10*.*  Recent  experiments 
by  Katz  et  al.  (1975)  performed  in  a diffusion  cloud 
chamber  show  that  the  nucleation  of  benzene  above  the 
triple  point  is  well  described  by  the  classical  theory. 

In  fact,  any  experimental  result  exhibiting  a 
"high”  nucleation  rate  in  such  comparisons  may  be  sub- 
ject to  the  additional  criticism  that  vapor  impurities 
have  resulted  in  binary  nucleation.  Such  nucleation 
yields  a lower  critical  supersaturation.  In  fact  the 
highest  values  of  J (or  D are  seen  in  Table  6.3  for 
chloroform  in  lines  6 to  11  based  on  nozzle  and  shock 
tube  experiments  performed  at  MIT  and  Yale  respectively. 
The  highest  rate  value  of  nucleation  T = 10**,  was  noted 
in  a shock  tube  (line  11)  whose  cooling  rate  was  lower 
then  that  of  the  nozzle  experiments  lines  6 to  10.  The 

liquid  chloroform  employed  in  the  experiments  at  our 
laboratory  (Barschdorff  1972b)  contained  0.75%  ethanol 
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as  a stabilizer.  Apparently  the  same  grade  (Russell 
1971)  was  used  previously  in  the  HIT  nozzle  experiments 
(Dawson  et  al.  1969)  seen  recomputed  in  Table  6.3, 
lines  6 to  10.  Condensation  appears  roughly  at  the 
triple  point  so  that  the  property  problem  does  not 
enter  in  this  case.  In  view  of  the  prior 
experience  with  binary  nucleation,  the  chloroform  data 
with  ethanol  present  in  the  mixture  cannot  be  accepted 

as  examples  of  homogeneous  nucleation.  Moreover,  recent 
experiments  by  Katz  (1975)  with  chloroform  using  a dif- 
fusion cloud  chamber  showed  nucleation  rates  close  to 
that  given  by  the  classical  theory.  Thi«  result  was 
obtained  with  both  the  sample  containing  U.7St  ethanol 
and  a purified  sample  believed  to  be  free  of  ethanol. 

The  data  on  chloroform  should  therefore  be  considered 
inconclusive  at  present. 

In  sum,  the  pure  vapor  expansions  tabulated  yield 
results  roughly  compatible*with  the  classical  theory 

of  nucleation  with  "low"  values  of  F.  Without  further 
defining  the  meaning  of  ”low"  in  view  of  all  the  uncer- 
tainties facing  us,  the  results  in  Table  6.3  tell  us 
that  the  substances  benzene,  ammonia.  Freon  11,  and 
carbon  tetrachloride  yield  high  rates  of  nucleation. 

Careful  additional  experiments  ruling  out  all  possibil- 
ities of  binary  nucleation  and  the  study  of  other  sub- 
stances will  be  required  to  further  verify  the  varying 
nucleation  rates  given.  It  is  interesting  to  note  that 
no  results  on  any  substance  in  either  Wilson  cloud 

* We  recall  that  employing  e.g.  Volmer's  equation  (Volmer 
1969)  which  is  essentially  identical  to  those  of  other 
early  investigators  and  viewing  T as  an  adjustment  factor 
as  in  Eq.  (6.1)  the  classical  theory  yields  T = 1. 

■ 1 


chambers  or  diffusion  cloud  chambers  have  yet  exhibited 
the  behavior  of  the  materials  cited  giving  "high**  values 
of  T.  It  is  thus  truly  remarkable  that  a theory  fraught 
with  such  uncertainties  provides  generally  consistent 
answers  on  nucleation. 


1 
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LIST  OF  SYMBOLS 

A Nozzle  cross-sectional  area, 

a Speed  of  sound,  Eq.  (2.10). 

b Eq.  (5.29). 

c Wave  speed,  Eq . C 2 . 37 ) t also  specific  heat,  Eq.  (5.10). 

cp  Specific  heat  at  constant  pressure. 

cy  Specific  heat  at  constant  volume. 

D,  De  Impingement  and  evaporation  rates,  respectively, 

Eqs . (5.33),  (5.34). 

d Lq.  (5.30). 

e Internal  energy. 

f(r)  Size  distribution  function. 

g Condensate  mass  fraction. 

h Enthalpy;  also  nozzle  height,  Eqs.  (2.33),  (2.34);  also 

Planck's  constant,  Eq.  (5.26). 

I Moment  of  inertia. 

J Nucleation  rate. 

K [=  (Y*1)/2(y-1)]. 

k Boltzmann's  constant. 

Kn  Knudsen  number,  Eq.  (2.3). 

L Latent  heat  of  condensation. 

I Characteristic  length. 

M Mach  number , Eq . (2.2). 

m Molecular  mass  (=  p/N^);or  index  of  reflection. 

N Number  density  of  droplets, 

n Number  of  molecules  in  a cluster. 
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Na  Avogadro's  number. 

p Pressure. 

Pr  Prandtl  number,  Eq.  0.13). 

0 Mass  flow  rate  through  a nozzle,  Eq.  (3.13). 

q Heat  addition  to  unit  mass  of  the  flow. 

qr€p  Replacement  factor,  Eq.  (S.76). 

R Molar  gas  constant t also  radius  of  curvature  of  nozzle 

contour,  Eqs.  (2.33),  (2.3«). 

r Droplet  radius. 

Re  Reynolds  number,  Eq.  (2.1). 

s Entropy ? or  saturation  ratio  (*  py/pw). 

T Temperature. 

t Time. 

u Velocity. 

v Molecular  volume. 

x Spatial  coordinate. 

a Mass  accommodation  coefficient. 

o'  Effective  condensation  coefficient,  Eq.  (S.35). 

T Eqs.  (5.25)  and  (6.1). 

Y Ratio  of  specific  heats,  Eq.  (2.9)  (=  Cp/cy). 

n Viscosity. 

k Thermal  conductivity. 

A Mean  free  path?  also  wavelength, 

y Molecular  weight, 

v kinematic  viscosity  (=  n/p)- 
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£ Thermal  accommodation  coefficient. 

0 Density. 

p'c  Density  of  condensate  referred  to  the  volume  of  the 

gas  phase. 

0 Surface  tension. 

t Relaxation  time. 

♦ Relative  humidity  * (pv/p^)  * 100%. 

^ Progress  variable. 

Y Eq.  (3.3). 

u>  Mass  fraction. 

Subsripts  and  superscripts, 
ad  adiabatic, 

c condensate, 

cf  classical, 

crit  Refers  to  critical  point, 

d Refers  to  droplet, 

e Equilibrium, 

exp  Experimental, 

f Frozen  flow. 

1 Inert  carrier  gas. 

iso  Isothermal, 

noniso  Nonisothermal . 

r Refersto  droplets  with  radius  r. 

s Shock  $ also  for  saturation. 
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th  Theoretical, 

v Vapor. 

0 Stagnation  conditions;  aiso  for  conditions  in  the 

supply. 

01  Supply  conditions. 

1,2, 3, 4 Refersto  shock  tube  flow  regions.  Fig.  2.2. 

• Refers  to  bulk  properties. 

* Sonic  conditions  at  nozzle  throat; or  for  critical 

clusters . 
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APPENDIX  1 

Compressible  Flow  Table.  Perfect  Gas,  y = 7/5  = 1.400. 


M 

P/P0 

P/PQ 

T/T 

o 

A/A* 

u/a* 

o 

o 

1.0000 

1.0000 

1.0000 

00 

0.0000 

CM 

• 

o 

0.9725 

0.9803 

0.9921 

2.964 

0.2182 

0.4 

0.8956 

0.9243 

0.9690 

1.590 

0.4313 

0.6 

0.7840 

0.8405 

0.9328 

1.188 

0.6348 

00 

• 

o 

0.6560 

0.7400 

0.8865 

1.038 

0.8251 

1.0 

0.5283 

0.6339 

0.8333 

1.000 

1.0000 

1.2 

0.4124 

0.5311 

0.7764 

1.030 

1.1583 

1.4 

0.3142 

0.4374 

0.7184 

1.115 

1.2999 

1.6 

0.2353 

0.3557 

0.6614 

1.250 

1.4254 

00 

• 

rH 

0.1740 

0.2868 

0.6068 

1.439 

1.5360 

O 

• 

CM 

0.1278 

0.2300 

0.5556 

1.688 

1.6330 

2.2 

0.09352 

0.1841 

0.5081 

2.005 

1.7179 

2.4 

0.06840 

0.1472 

0.4647 

2.403 

1.7922 

2.6 

0.05012 

0.1179 

0.4252 

2.896 

1.8571 

00 

CM 

0.03685 

0.09463 

0.3894 

3.500 

1.9140 

3.0 

0.02722 

0.07623 

0.3571 

4.235 

1.9640 

3.5 

0.01311 

0.04523 

0.2899 

6.790 

2.0642 

4.0 

0.006586 

0.02766 

0.2381 

10.72 

2.1381 

4.5 

0.003455 

0.01745 

0.1980 

16.56 

2.1936 

5.0 

0.001890 

0.01134 

0.1667 

25.00 

2.2361 

10.0 

0.0000236 

0.000495 

0.04762 

535.9 

2.3905 

APPENDIX  2 

Condensation  Studies  in  Compressible  Flow  Systems. 

A chronological  list  of  publications  on  homo- 
geneous nucleation  in  flow  systems  is  assem- 
bled. These  papers  contain  experimental  data  and/or 
they  are  based  on  the  authors ' own  or  other  experimental 
results.  Investigations  motivated  by  the  under- 
standing of  the  physics  of  homogeneous  nucleation  as 
well  as  those  in  which  technological  problems  inspired 
the  studies  are  both  listed.  No  claim  is  made  for 
completeness.  In  fact  owing  to  the  uncertainties 
attached  to  many  "air"  condensation  studies  in  hyper- 
sonic wind  tunnels  (i.e.  wind  tunnels  for  M>>1),  few 
experiments  on  "air"  condensation  are  included.  More- 
over the  many  low  speed  condensation  studies  on  two- 
phase  equilibrium  flows,  studies  of  purely  diabatic 
flow  and  of  unsteady  nozzle  flow  with  condensation  are 
omitted.  Finally,  the  field  of  clustering  in  molecular 
beams  (see  Section  4c)  is  completely  omitted;  it  was 
recently  covered  extensively  by  Hagena  (1974),  and  it 
falls  outside  the  scope  of  the  current  article. 
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TABLE  CAPTIONS 

4.1  Legend  for  Fig.  4.4.  Collection  of  onset  pressures  and 
temperatures  for  condensation  of  water  vapor  in  steady 
supersonic  nozzle  flow. 

6.1  Number  of  water  and  ethanol  molecules  respectively  in 
spherical  drops  of  a given  radius  computed  with  Eqs. 
(5.22)  and  (5.23).  Density  for  H20  and  CH3-CH2-OH 
condensate,  1 and  0.86  g/cm3 , respectively,  T = 210  K. 

6.2  Comparison  of  condensation  results  in  nozzle  and  shock 
tube  flow  with  nucleation  theory.  Data  given  pertain 

to  the  state  at  the  onset  of  condensation  and  the  supply 
of  pure  vapors  (steam  and  N2 ) . (For  theoretical 
expressions  used  note  Sections  5 and  6.) 

6.3  Comparison  of  condensation  results  of  vapors  in  inert 
carrier  gases  in  steady  nozzle  flow  and  shock  tube 
experiments  with  nucleation  theory.  Data  given  pertain 
to  the  state  at  the  onset  of  condensation  and  the  supply. 
(For  theoretical  expressions  used  note  Sections  5 and  6.) 
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FIGURE  CAPTIONS 


1.1  Schlieren  picture  of  a "condensation  shock"  in  a super- 
sonic nozzle  operated  with  moist  air.  This  picture 

taken  by  A.  Busemann  and  shown  at  the  Volta  Congress  in  Rome 
by  Prandtl  (1935)  is  believed  to  be  the  first  published  photo- 
graphic record  of  condensation  in  supersonic  flow.  The 
"Mach  lines"  emanating  from  the  walls  are  caused  by  the 
purposely  roughened  nozzle  surface. 

1.2  Schematic  p-T  diagram  with  coexistence  line  and  isentrope. 

The  adiabatic  supercooling,  AT^ , and  the  critical  super- 
saturation, P^/Poo  are  defined. 

1.3  Adiabatic  supercooling  of  water  vapor  in  air  as  a function 
of  relative  humidity  in  the  atmospheric  supply  of  a 
(40x40 )cm2  supersonic  wind  tunnel  in  steady  flow.  Data 
obtained  from  Schlieren  photographs  (Eber  and  Gruenewald, 
1941/42)  with  a M = 1.86  nozzle  at  a cooling  rate  of 

= 105  c/sec. 

1.4  Adiabatic  supercooling  of  water  vapor  in  air  as  a function 

of  cooling  rate  in  a steady  Prandtl-Meyer  expansion 

occurring  in  supersonic  flow  around  a sharp  40°  corner, 

for  41  < $ < 69%  (Smith  1971).  The  solid  line  shows  an 

o 

interpolation  of  the  results  with  the  limit  for 
dT/dt  = 0. 

2.1  Area  ratio  and  dimensionless  flow  properties  as  a function 
of  Mach  number  in  steady,  inviscid  nozzle  flow  with 

y = Cp/cv  = 7/5  = 1.4  (e.g.  air). 

2.2  Schematic  diagram  of  the  ideal  operation  of  a shock  tube. 

At  top  a time-distance  plot.  At  time,  t^ , the  instan- 
taneous pressure  and  temperature  is  shown  below  with  the 
corresponding  flow  features  in  the  tube  given  at  the 
bottom. 


1 


129 

4.1  Experimental  pressure  distribution  in  two  nozzles  as 

a function  of  distance  in  steady  nozzle  flow  from  an 
atmospheric  reservoir  or  supply  with  and  without  con- 
densation (Wegener  and  Pouring  1964).  (a)  4>o  = 29%, 

dT/dt  = - 0.7  C/ysec.  (b)  4>q  = 34%,  dT/dt  = - 2.8  C/ysec. 

The  symbol  indicates  the  results  of  a new  complete  solu- 
tion of  the  equations  of  motion  together  with  nucle- 
ation  rate  and  droplet  growth  with  log^T  = 4.7  and 
a = 0.6  (see  Sections  5 and  6). 

4.2  Flow  parameters  computed  from  the  diabatic  equations 

of  motion  through  a heat  addition  zone  caused  by  con- 
densation of  water  vapor  in  air  in  a steady  nozzle  flow  in  a 
(18x18)  cm2  nozzle  for  M = 2.92  (Wegener  1954).  Pressure 
distributions  were  measured  for  dry  air  and  an  atmo- 
spheric supply  with  4>Q  = 53%.  and  po,  are  supply 

temperatures  and  pressure  with  Tq  and  pQ  being  local  values, 

respectively.  Experimental  onset  of  condensation  is 
observed  at  x = 1.6  cm. 

4.3  Mean  water  cluster  size  and  concentration  as  a function 
of  initial  specific  humidity  (H2O  mass  fraction)computed 
from  light  scattering  measurements  assuming  a delta- 
function  size  distribution  as  computed  from  Eqs.  (4.1  ) and 
(4.2)  (Stein  1969a).  Steady  supersonic  flow,  atmospheric  supply 
of  moist  air,  thermodynamic  equilibrium  at  the  obser- 
vation station  downstream  from  the  condensation  zone. 

Pressure  and  temperature  at  the  onset  of  water  conden- 
sation (or  critical  supersaturation)  in  expansions  of 
steam  and  moist  air  in  steady  supersonic  nozzle  flow. 

Different  initial  states  and  cooling  rates  are  shown  as 
identified  in  the  references  listed  in  Table  4.1  giving 
the  symbols  of  this  figure. 

Schematic  of  symbols  and  states  used  in  droplet  growth 
calculations . 
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6.1  Observed  (shaded  area)  and  computed  (solid  line)  mass 
fraction  of  liquid  ethanol  in  air  in  steady  supersonic 
nozzle  flow  (Wegener  et  al.  1972).  The  shaded  area 
indicates  the  experimental  uncertainty.  Nucleation  rate, 
number  density  of  clusters  and  saturation  ratio  as 
computed  for  this  experiment  with  an  atmospheric 

supply  and  <t)Q  = 0.005  . Note  distance  and  area  scales. 

Time  scale  computed  with  Eq.  (2.32).  In  the  calcula- 
tions with  nucleation  rate  and  growth,  T = 10s,  a = 0.7 
and  a = a , see  Section  5. 

6.2  Critical  radius  and  selected  droplet  radii  in  steady 
nozzle  flow  for  liquid  ethanol  as  computed  from  Volmer's 
Eq.  (5.20)  with  a/a oo  = 0.725  for  the  experimental 
environment  described  in  Fig.  6.1.  Normalized  drop 
size  distribution  at  the  nozzle  exit  is  shown  on  the 
right . 

6.3  Experimentally  observed  critical  supersaturation  as  a 
function  of  temperature  for  steam  condensation  in  steady 
supersonic  nozzle  flow  and  a shock  tube.  The  solid  line 
shown  is  derived  from  all  supersonic  nozzle  steam  experi- 
ments listed  in  Table  4.1.  Temperature  made  dimension- 
less by  Tcrit  = 647  K. 

6.4  Ratio  of  observed  and  computed  nucleation  rates  (non-iso- 
thermal  theory,  Eqs.(5.20)  and  (5.28)  at  the  onset  of  condensa- 
tion as  a function  of  initial  specific  humidity  obtained 

in  the  expansion  fan  of  a shock  tube  (Barschdorff  1975). 

Note  different  carrier  gases:  O air,  □ helium,  O argon, 

£ none  (steam) . 

6.5  Ratio  of  observed  and  computed  nucleation  rate  at  the 
onset  of  condensation  of  steam  in  supersonic  nozzle  flow 
(Barschdorff  et  al . 1972).  Results  from  non-isothermal 
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nucleation  theory,  Eqs.(5.20)  and  (5.28) , and  isothermal 

theory  Eq.  (5.20) , are  given  as  a function  of  conden- 

sation coefficient  for  a single  experimental  result  for 
P/Poo  = 8,1  (Barschdorff  1971),  cf.  line  3,  Table  6.2. 

6.6  Critical  supersaturation  as  a function  of  temperature  of 
water  vapor  present  in  small  concentrations  in  carrier 
gases.  Temperature  function  according  to  Eq.  (6.4). 

Results  from  cloud  chambers,  nozzles  and  shock  tubes  are 
shown. 

6.7  Critical  supersaturation  of  steam  from  nozzle  and  shock 
tube  experiments  plotted  as  in  Fig.  6.6. 

6.8  Matched  values  of  the  surface  tension  and  the  replacement 

factor,  Eq.  (5.20)  compatible  with  experiments  on  the 

onset  of  condensation  of  ethanol  in  air  (w  = 0.005)  for 

o 

the  experiment  described  in  Fig.  6.1  (Wegener  et  al . 1972). 
Error  bar  indicates  the  uncertainty  in  the  measured  nucle- 
ation rate. 

6.9  Numerical  values  for  the  terms  entering  the  nucleation  rate 
Eqs.  (5.20,  (5.25),  and  (5.26)  in  steady  nozzle  flow  of  ethanol 
(w  = 0.005)  in  air  of  the  experiment  described  in  Fig. 

6.1  (Wegener  et  al.  1972).  Independent  calculations  of 
the  replacement  factor  in  the  expression  of  Eq.  (5.26) 
for  T,  following  the  models  of  Lothe  and  Pound  (1966) 
and  Dunning  (1965)  respectively.  Computed  experimental 
value  of  T,  based  on  Eq.  (6.1),  see  also  Table  6.3. 


Agreement  of  Te  and  implies  that  the  classical  nucle- 
ation rate,  e.g.  Eq.  (5.20)  gives  results  remarkably  close 
to  the  experiment  at  the  onset  of  condensation. 


TABLE  4.1 
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Symbol  Author  (year)  Cooling  Rate  (C/ysec) 


<? 

Binnie  6 Green  (1943) 

1 

A 

Pouring  (1963) 

0.7 

and  3 

□ 

Gyarmathy  6 Meyer  (1965) 

1 

and  0 . 3 

X 

Stein  (1967) 

2 

b 

Deych,  Stepanchuk  6 Saltanov 

(1969) 

2 

▼ 

Roberts  (1969) 

1 

• 

Barschdorff  (1971)  and 
Barschdorff  6 Ludwieg  (1972) 

0.6 

+ 

Stein  6 Moses  (1972) 

1 

